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INTRODUCTION 



I . 

Recent developments in control theory and practice 
coupled with digital electronics advances have produced 
increased interest in digital power plant modeling. Although 
conventional marine boiler modeling is not a new concept, it 
is still an area of increased interest due to the relative 
lack of analytical data. This deficiency is due primarily 
to the extreme environment that must be endured by sampling 
equipment. Advances in data collecting and processing equip- 
ment have made boiler data collection practically feasible, 
and in fact, nuclear steam generator development depends 
heavily on actual data. However, it is not economically 
feasible for a single manufacturer to retrofit modern data 
collection systems on a standard marine boiler. The majority 
of marine boiler models available are not attractive from 
the control engineer's viewpoint because of one or more of 
the following: 

a. The model is over-simplified to the point of being a 
"teapot" model. 

b. The model is completely developed and written in 
Laplace transforms and the state space equations are 
too difficult to extract. 

c. The model is over-complicated to the point of being 
computationally inefficient. 
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d. The causes and modeling of "shrink and swell" 
phenomena are ignored. 

This model attempts to compromise between simplification 
and complication while including a theory and model for shrink 
and swell. The model follows the path specified by Fini [1] 
restated below: 

a. A general D-type marine boiler model is prepared in 
Continuous System Modeling Program (CSMP) language. 

b. A general FORTRAN program is used to prepare the 
initial conditions needed by the CSMP model. 

c. The initial condition program depends only on data 
easily obtained from the manufacturer's technical 
manual and engineering handbooks. 



, 
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II. MODEL CONSIDERATIONS 



A. BOILER TYPE 

The D-type boiler was chosen for modeling because of 
its comparative standardization among manufacturers coupled 
with its popularity in naval ship power plants. This 
particular boiler is one installed in the LHA-1 class of 
United States naval ships. Since D-type boilers operate 
under a wide range of geometry and output conditions the 
initial condition and model programs are designed to operate 
over similar ranges. 

B. MODEL NOMENCLATURE 

The model nomenclature differs from reference I 1 ] because 
of a need to clarify the model equations and afford readability 
to the actual computer programs. The following objectives 
were considered when developing the nomenclature. 

a. The equations must be easily read, necessitating 
minimum variable-definition referrals by the reader. 

b. The notation used in the model development must be the 
same as that used in the computer programs. 

The resulting notation combines a. and b. and consists 
of frequently used notations for variables followed by a 
subscript to indicate state point location and/or time using 
Figure 1. 

Table 1, a tabulation of general rules and examples, 
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should suffice to allow the reader to follow the model 
development . 

C. EQUATION DEVELOPMENT FORMAT 

The dynamic model and initial condition programs are 
developed simultaneously with corresponding equations 
written in appropriate proximity to each other. A notation 
of [IC] following an equation indicates that it is an initial 
condition equation, while a notation of [ME] indicates a 
model equation. In some instances the equation applies to 
both programs, an indication of [IC, ME] is used to indicate 
this . 

In as many instances as possible equations are developed 
directly from elementary heat transfer and fluid dynamic 
principles. This stems from a desire to encourage more re- 
search in this particular field and to foster a better under- 
standing of marine boiler dynamics among students of marine 
engineering. The model was developed to be read and used by 
students with a minimum of additional research. 

D. ASSUMPTIONS 

Various assumptions are made to facilitate either model 
simplification or program tractabil ity . The latter are not 
always desirable and further research or different programming 
techniques may obviate the need for them. The following 
major assumptions are made: 

1. The energy transferred to the generating tubes is 
distributed uniformly, i.e. "even heating" exists. 
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2. There is a thorough mixing of the fluid in the steam 
drum. 

3. Steam generation occurs only in the boiling section 
of the screen and main bank risers. 

4. All generating tubes and downcomers in a particular 
circuit have the average geometry. 

Uniform energy transfer is assumed to facilitate model 
simplification and is justifiable since the generating 
banks are designed to optimize uniform heat transfer. Mixing 
of the fluid in the drum is a result of the turbulence in the 
drum. The third assumption results from a lack of reliable 
experimental data and the knowledge that D-type boilers are 
natural circulation boilers relying on a density difference 
in the circulation loop to promote flow. Steam production 
outside of the generating banks would decrease the density 
difference, thus inhibiting natural circulation. The geometry 
assumption is necessitated by model simplification. 






12 



III. MODEL DEVELOPMENT 



A . GENERAL 

The boiler is developed from control volumes for which 
the appropriate equations for energy transfer and storage 
rates are written. The initial condition equations are 
generated by setting the rate (d/dt) terms equal to zero and 
solving for the desired variables. 

The model is developed with three different elements: 

1. Furnace-side heat transfer 

2. Water-side heat transfer 

3. Water-side circulation 

B. FURNACE-SIDE HEAT TRANSFER 

The gas flow path can be followed on Figure 1 where it is 
represented by the dashed line. The fuel and air mixture 
enters the furnace at state point Q where radiation heat 
transfer occurs to the screen riser tube metal. The flue gas 
then leaves the furnace and exchanges heat via convection with 
the superheater (state point R-S), main bank riser tubes 
Cstate points S-T), and economizer (state points T-U) in that 
order . 

At state point Q the fuel-to-air ratio and fuel and air 
temperatures are relatively constant. This allows the 
sensible heat of the fuel and air to be lumped. The mass flow 
rate into the furnace is the sum of the air flow and fuel flow. 
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M<3Q = MFQQ -hjrviAqQ [>&][ 

The initial values of fuel and air flow at a specified 
operating point are known. Therefore, the corresponding 
initial condition equation is: 

The mass flow rates through the remaining sections of the 
boiler are the same. 



Mq<a = W(iUL= mss = rwrr- rnu.u. £ M e] 

= M*1^ = MSS^ = m7 Tfi = MuUf f 

The total heat supplied to the boiler can be written as 
the product of the mass flow rate of fuel and the "lumped" 
fuel heating value. 



G[G! = MFQ3 * FtfV 






At a specific operating point the heat released to the furnace 
is available from the technical manual as is the furnace 
volume. The total heat input at steady state is the product 
of these two. 

* fusvol. ^re] 



This allows the computation of the fuel heating value. 



Ftt\/ - ^xcTj 
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1 . 



Furnace to Screen Risers 



The heat transfer to the furnace screen tubes occurs 
primarily through radiation, i.e. 

Gj - CT*A ( lgg, -Tw) 

which is the Stef an-Boltzmann law^where a is the Stefan- 
Boltzmann constant defined as: 

- a.?sv7 x )0 ~ n am 

Therefore the heat transfer to the furnace screen can be 
written 

Ql = 

~ (t W + * *. q.0 ) £meJ 

where 

Si h 

At a given operating point, technical manual values are 
specified for the furnace heat absorption rate and the area 
of the radiant heat absorbing surface of the furnace screen. 
The product of these two yields the steady state heat transfer 
rate which in turn can be used to compute SIGMAA. 



This equation assumes that all the radiation from the combus- 
tion flame strikes the screen tubes, and both the flame and 
the tubes behave as black surfaces. 
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QL0 = G)Fuj?A 6 4*. Al^AO 






Si^MAA = q $4 / C(trr 0 -i-4k0.0) * * 4.0 - 
(TWfi + 460.0) * * 4.0) 

An energy balance equation for the mass of flue gases in the 
furnace region can be written in the form 

-d ( yy\ dp T) - 4d*al ^upplt'cd — energy 

cH" "lrdnSTCrrej "ho Screen 4-uJoe.s — 

Srieropj o*f l©3*AV)(j "furnace- 
rao.; 0 r\ 

or ^ 

d. (wAssqe. * cqfi. * tw?.) = qq -qi- m r/z 
dt * cqe. -sfc (nse.-TArvi6) 



Note: TRR-TAMB is the absolute temperature. At steady state 
the rate terms are equal to zero. Values of total heat 
supplied to the boiler, heat transfer rate to the screen 
tubes, total flue gas flow rate, and flue gas temperature 
are available from previous equations. This allows the specific 
heat of the flue gas in the furnace region to be computed. 

- (qqj i - <310 ) / % (tqRj 0 -TP/*fi)) \x<f] 

During a transient, the density of the flue gas in the furnace 
and the specific heat of the furnace flue gas can be regarded 
as constants near a steady state operating point. The dynamic 
equation can be written as 
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OTRfS. - (qQ -Ql- MfcJZ. ■*. CQfc. ^ (l~RR. ~~ TAN\&}) / 

(r'vAssqie. - jk. cqi^,) \ me ] 

where 

NASSQR. - fcFLuE. ^ Fui^VOL Yxc] 



Thus , 

Tfilt 





Orae -t- 



or in CSMP-III language 



~ ifora^u^TRR^ > dtrjz ^ 




The value for TRRO is available from the technical manual. 

2 . Flue Gas to Superheater 

For the remaining tube banks heat transfer occurs 
primarily through convection. The energy given up by the 
flue gas in flowing over the tube surfaces is represented 
by the general equation 



energy = mass flow * specific 

given up rate of heat of 

flue gas flue gas 

* change in temperature of the gas 



i . e . 



°1 



v'x (L r 



M 
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At steady state, the energy transferred to the superheater 
tubes from the gas must equal the energy transferred from 
the tubes to the steam. 

The steady state values for superheater inlet and outlet 
steam temperature and pressure are specified by the technical 
manual along with the superheated steam mass flow rate. This 
allows the direct computation of Q40 and Q30. 

(l4UNJ^ - HMM0) [icj 

HMMO is the enthalpy of saturated steam corresponding to drum 
pressure, computed with curve fitted equations [2]. HNNO is 
available directly from the technical manual. The specific 
heat of the flue gas in the superheater region can now be 
computed . 

CRS ^ (TRRJ0 -TSS0V) 

Thus the dynamic equation for heat transferred from the flue 
gas to the tube metal is written: 

Q3 - ^ s (trr-Tss) [>\eQ 

The heat transferred to the tubes via convection may also be 
represented by the Grimson correlation £ 3 J which states that 
for cross flow over tubes; 
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" Pr V3 



^cL - 

H 




where 

d 

<L 

Pr 



Ve.loc.4y 

"Via. loo cti dn^si^r^ 
C>oia^Ow'i~ 

Pra^cl'H n tcrr. b&r' 



Since h can also be written 

K - c| /\ST 



then 






l\<L 




n 



also 

"UeO — -A. 

•?* 

leading to the general equation 




Pr 1/3 AT 



Pr 1 ' 3 



Around a particular steady state operating point the physical 
properties can be regarded as a constant and the physical 
dimensions are constant, i.e. 



AC, 




Pr 1/3 = Constant 



which implies 



Cj - fcons+anO n 



For flow over tube banks with differing rows of tubes and 
differing tube geometries the constant (n) varies little and 
the average is approximately .6 [3]. Therefore, in the case 
of the superheater and main bank riser tubes, the correlation 
is used in the form 



Since the steady state values for Q30, MRRO, TRRO , TSSO, and 
TWWO are given by the technical manual the constant KRS can 
be calculated. 



KftS •=. Q 3 ?z5 / ^ .(a * (TRS0 -TVJWjtf)) fee] 



C| - Cr^-^AT 



or for the superheater 



Q3 - KRS * MRR. * .6 (tR-S — TWW) 



where TRS is the average flue gas temperature 
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To provide the dynamic solution of superheater flue 
gas outlet temperature, the two transient equations for the 
flue gas to superheater tube metal heat transfer are equated 
resulting in 

"rfts - rww - ^ crs /krs 

Tfcfc - TSS 

where 

Tfcs - (rgjz. tss} /z.pS 

Solving the above for TSS gives 

Tss - (trr, £™ e ] 

where 

PH-n ~ tK fV\R*L^ * ^ C^s/kRS Yme] 1 

3 . Flue Gas to Main Bank Risers 

The equations for the flue gas to tube metal energy 
transfer in the main bank are derived similarly. For the 
equation involving flue gas specific heat; 

~ rviss cst ^ (tss — t rr) 

The Grimson correlation equation for the main bank risers is 

- KST mss ik ik 0 .(q ^ (tsT - TvV) 
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where 



1 ST = (TSS +TTr) /Z.0 

By equating these two equations as was done with the super- 
heater the flue gas temperature leaving the main bank can 
be written 

ttt = (tss % (pHrz-l)+2.^Tw)/(PHr^+l) LmeJ 

where 

PHIZ ~ Z,* hSS * 0.4 ^ CST /k^T 

The steady state heat transfer is computed differently than 
that for the superheater. Since the steady state heat trans 
fer to the screen bank has been computed, the energy 
transferred to the main bank is the difference between the 
total energy transferred and the desuperheater and screen 
bank energy transfer. At steady state the mass flow rate of 
steam out of the boiler is equal to the mass flow rate of 
feed in the boiler and the total energy transferred is the 
mass flow rate multiplied by the enthalpy change between the 
outgoing steam and incoming feedwater. Thus, the steady 
state energy transfer equation for the main bank is 

G{<b$ = MmMj z$ t rc 0 

Q990 is computed in the water side heat transfer section. 

Since Q50 = Q60 at steady state the heat transfer 
coefficient KST and the specific heat CST can be computed. 

KST - / (mss0 ^ 0.6 ^ (tst0 -tyy 0)) £rcl 

csr -350/ (rv\ss0*. (rss0-rrr0)) \rc] 



22 



4 . Flue Gas to Economizer 

The economizer tubes differ from the generating and 
superheater tubes in that they are finned. A reasonable 
finned tube heat transfer correlation was developed by 
Weirmann, et al . [4] 

j = V, Pr 2/3 

Cp G.r vxa* 

G s maximum mass flow rate 
max 

j = correlation constant 



Since 




specific heat of flue gas 




the above 



correlation can be written 

C| ~ j A AT* CjQ G. max 
P 



G is a function of mass flow rate into the tube bank and 
max 

tube bank geometry. Similar to the case of the Grimson coeffi- 
cient used for the superheater and main bank tubes the physical 
and geometric properties are considered constant leading to 
the use of the correlation in the form 

- Cm i) 
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Therefore, for the economizer the equation is 

Q7 ~ KXU. tK iv\tt ^ (ttul -tkx) 

where 

TTUL - TTT -h TU.U. 

Z.0 

Paralleling the superheater and main bank, the heat transfer 
may also be written 

Q7 = MTT CTU. ^ (ttt-Tu. [meT) 

The steady state equation development follows that of the 
superheater and main bank sections and is not repeated here. 
The equations are: 

- MAA0 * (HSeszS-HAAjzO 
67jz( = Qsy tiaj 

ctu. ^ Q,lfi/(mr<p -Tu.iqzS)) [rq] 

Kru. = G70/ (MTTjZS *(TTU.0 - £ r c.] 

TTUjz5 = (ttt^ -t- Tuufi) /z.jzJ ^xe] 

Knowing the values for KTU and CTU developed by the initial 
condition program allows the solution of the flue gas outlet 
temperature in the transient analysis. Thus, 

ime] tuu - (ttf #. (pHr3 %.<f > * tm)/ 

(PKX3 -t-l.0) 

Ime: PH-I3 =2.{zS^ CTU./krL<_ 
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C. WATER SIDE HEAT TRANSFER 



1. Risers 



The heat supplied to the riser banks is used both in 



the boiling and nonboiling length of the tube; however the 
assumption is made here that the heat transfer is all under 
fully developed flow boiling conditions. Since the nonboiling 
length of the risers is relatively small, little error is 
produced. For fully developed nucleate flow boiling, Levy 
[5] suggests the equation 



By lumping the constant area and denominator the riser equations 
can be written in the forms 



- kv PsArT * ^ (4. < 75 / 3 . 0) ik (rW- tsjvt)*. -K 3.0 ImeJ 



Q(q - KV * fs/vr 95/3 .^)^ (ttv -t svr)^ *3.0 



for the screen and main bank respectively. A simple energy 
balance on the tube metal yields the metal temperatures, i.e., 





3 



Therefore, the two riser equations are 



and 
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The equations for screen and main bank metal temperatures are 



Otv v ■=• (q t -Qz)/ ( m ftss v epv) ej 

TV\J - iMrSfiL , OT/i/) [Vl£~] 

Orvv - (qr-qfc)/ ( i^assy ^ CPy 1 ) [we] 

tvv =• ikjtssu-(tyy{zS ^oryv) 



The heat transfer coefficient KY can be computed 
directly from tube geometrical data available in the 
technical manual ; 

KY - AieeAMe/^ga-E. jzS<b Exc] 

Since Q10 = Q20 at steady state, the dynamic equation for 
Q20 can be solved in a steady state situation for KV, 

tic) 10/-= qi^PSATVf* ( 4 / 3 )* (rW0-TSffl)**3.?5) 

Q60 was computed in Section III-B.3. This allows 
the calculation of the initial main bank riser tube metal 
temperature , 

[ ic] TYitf = ( 'Q(o 0 /(Mk PSfirT*.* (4/3))) ^ 6-/3) -t-TSPtT 

In the case of the screen risers, the initial tube 
metal temperature is available from the technical manual and 
the initial heat transfer rate has been previously calculated 
This allows the calculation of the screen riser steamside 
heat transfer coefficient, 

tic] fw = qijz 5 /(PSAT * * (m.(zS/3.0)^ 3.0) 
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2 . Desuperheater 



The Dittus-Boelter correlation [6] is used to compute 
the heat transfer rate from the steam to the desuperheater 
tubes. This correlation in its basic format is 

NJu. - 

where 

Nu => Nusselt number 
Re => Reynolds number 
Pr => Prandtl number 
n => .3 (for cooling fluid) 

When the appropriate variables and constants are substituted 
for the dimensionless numbers, the Dittus-Boelter correlation 
can be written 

0, - K.A AT 



where the constant K is defined as; 

^ ^ -Hienv'al Cottduuckv.ly ^ Wt- Vd^fec ansa 

-VuEo. 






( d y ? 

\*TT^: dU&nxei'er* ^ V\Scosi4y /K nr cE TuEesy 



* Pr* 



and AT is the log-mean-temperature difference (LMTD). 

Thus, for the desuperheater, the equation is 
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Q°[ - KfJ 9 |V\MfJ & ^ .? ^ LMTOfOP 



KNP can be computed directly. Thus, 



KlOP- ((.faz ^THCoiOfO)/oMP)^<(n.^/(Pr^DfOP 



7& VISGdkJ * KJTU80S))^K (Z i 
AfcEADS 



(=KAKJ eK* 0*3 

IxcQ 



The energy given up by the steam to the tubes may 
also be written in terms of steam specific heat, flow rate, 
and temperature difference. 

- d^P ^ -tap) 



Similarly since the heat transfer to the water in 
the drum is totally by convection; 

QSS = Kg *.(jgg -THH) t>\Ej 

For initial condition calculations THK is considered 
to be equal to TSAT. Therefore 

Kg = /(TZZ 0-TSKr) [i<J 



The desuperheater outlet temperature is solved by 
equating the two dynamic equations for Q9 and solving for 
the desired outlet temperature. 



tpp = (tkjm - Tggy^exp^p/ (e,tj p * <p.z ))) 

+TS - 5 Ime] 
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The instantaneous desuperheater tube metal tempera- 
ture is obtained by an energy balance on the tube metal 
similar to that for the riser tube metal. 



r*s ass K dT 

dt 





Solving the appropriate equation for dT/dt results in 

OTtt - tK CPs) \ja£J 

- IfOTSau DT&5- 



The specific heat of the steam in the desuperheater 
(CNP) may be calculated using the known steady state values 
for desuperheated steam flow rate and the desuperheater inlet 
and outlet temperatures and pressures. 

- MWW0 H PP0) 

ttJP - qS^>/(MtJiO9S^.(TWM^-TFP0)) [re} 



With the value of Q90 computed the initial log-mean- 
temperature difference can be calculated. 

LmtdiOP - q^/CMP 0.?) 

The log-mean-temperature difference is a function of 
steam temperature in, tube metal temperature, and steam 
temperature out; therefore the initial tube metal temperature 
may be calculated. 
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T - EKPOOS KTffrf)/ (\.<p -EXjpoDS^ [pi] 

exfODS- EXP((TKJkJ 5Z! ~TPP0 V L,v ' T OlOP) trcK 



3 . Economizer 

Paralleling the desuperheater development the Dittus- 
Boelter correlation is used to relate the heat transfer from 
the tube metal to the feedwater. Thus, 

tOu - .0ZZ R<?'*Pr n 

the only difference being that n is now .4 vice .3 because the 
fluid is now being heated. The appropriate constants are 
again lumped yielding 

Q? - tOC X MAA X X X LMT0A6 

An energy balance on the water gives the heat absorbed in two 
other forms ; 

■= K l AA ^ {h&S - HAa) 

Q? - x ca£ * ( t&q -taa) 



At the specified operating point the inlet and outlet condi- 
tions of the economizer as well as the mass flow rate are 
given, allowing the computation of the steady state heat 
transfer . 

§*<p = MAA$z( *. -HAA0) ire.] 
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The heat transfer coefficient KX is computed directly. 



kX= ((.^Z3 *THCon;a)/oA6)* (m.^/(pE^0A8tK 
VISCOA * NJTu6EC))^.3i< <f>. ? *. PRAA** ^-4 
*. AREA EC. £ x <0 



This permits the solution of the steady state log-mean-tem- 
perature difference. 

UwroAS - q SyS /(KX * rAAA^**.<zi.?) fr<CI 

Following the identical path delineated for the desuperheater, 
the tube metal initial temperature can be computed. 

TXX.JZ i = (tAA0-T6iS(Z$ * -EXfOFC.) [xcj 

EXPOEC = exp((t8s^-taa^ V LMTOAQ ltd] 



An energy balance on the economizer tube metal yields 



its instantaneous temperature. 

dtxx - (q?-G)7)/(mAssy. * c,Pa) [>e] 

TXX — t MTAitL. 5 OTXX. ) e J 



By equating the formulations for Q8 involving LMTD 
and specific heat, the economizer outlet temperature may' be 
calculated . 

T66= (taA -TXx)/(exP (kx/(cAS*.MAA#.* 0.z))) 
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The specific heat of the feedwater in the economizer 



is calculated at steady state using the computed value of 
steady state heat transfer and the given mass flow rates and 
inlet and outlet temperatures. 



4 . Superheater 

The development of the superheater equations follows 
that of the desuperheater and economizer and for that reason 
the development will not be repeated. The equations are 
listed below. 



The initial superheater log-mean-temperature difference 



is computed directly allowing the subsequent calculation of 
the tube-metal-to-steam heat transfer coefficient used in the 
Dittus-Boelter equation. 



Lr^TDmkJ - -TW0')/ ((TWW0 -Ti^M0) / 
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An energy balance on the tube metal yields the instantaneous 
tube metal temperature, the initial tube metal temperature 
being available from the technical manual. 

dtvJW - (q3-q 4) / (mASSw^ cpw) 

TWW - Orvjw) \L™ e I1 

The superheater outlet temperature can now be calculated. 

TfOlO - (rn^l -TWIa/) / (EXP (kW / (c,P\*J ^ 

^*^0)) VTWu/ ^jV\e] 

D. WATER-SIDE CIRCULATION 
1 . General 

The water/steam side circulation equations are by far 
the most difficult to visualize and codify. These equations 
must be accurate in order to predict phenomena such as shrink 
and swell while appropriate assumptions and simplifications 
must be made in order to make the equations tractable. 

As can be seen in Figure 1 , the feedwater enters the 
economizer at state point A, passes through the economizer 
to the steam drum where it becomes part of the water volume. 
The liquid leaves the drum via the downcomers, both main bank 
and screen at state points G and C respectively. The main 
bank downcomer delivers its liquid to the main bank risers 
via the water drum, and the main bank risers deliver the fluid 
to the steam separators. The flow through the screen risers 
is the same with the exception that there is no water drum in 
this circuit. The steam separators separate the majority of 
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the water from the steam leaving the drum at state point M 
and the majority of the steam from the water being discharged 
back into the drum liquid. The steam, with a negligible 
amount of water passes through the superheater via state 
point M-N . At the outlet of the superheater the steam leaving 
at path III is used for "main" steam; that leaving via path 
II travels to the desuperheater where it is cooled, then 
leaves via state point P and is used for auxiliary steam. 

The water leaving the separators with a small amount 
of entrained vapor is mixed with the incoming feed water and 
forms a "foaming" vapor/liquid mass in the bottom half of 
the drum. The liquid from this mass leaves via the downcomers 
and is circulated through the loop. 

2 . Downcomer Pressure Drop 

A closer look at the downcomer flow is necessary to 
justify assumptions that are made in the development of the 
pressure-drop equations. 

Circulation ratio is defined as the ratio of the total 
weight of steam liberated to the drum [7]. For a 600 pound 
marine boiler this circulation ratio is on the order of 20:1 
I 7 and 8] . This implies that for every 21 pounds of liquid 
flowing down the downcomers, 20 pounds of it has already 
traveled up the riser and is at saturation temperature. There- 
fore, assuming the downcomers are perfectly insulated, it is 
reasonable to assume that the downcomer liquid is at or very 
near saturation temperature. 
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The momentum equation for the downcomers may be 

written : 

pressure at top - pressure at bottom 
= frictional loss - gravitational head 
+ entrance loss + bend loss 
+ exit loss + inertia force 

The inertia force term may be considered negligible [9] . 

This is the quasistatic approximation which basically states 
that pressure waves move much more rapidly through the system 
than the important time constants. This is mathematically 
equivalent to the elimination of a large negative eigenvalue. 
This quasistatic approximation is only good for on-line 
transients and does not apply for extremely large discon- 
tinuities which would result in the boiler being taken out of 
operation, (e.g., a ruptured tube). Therefore the downcomer 
momentum equation can be written: 

pressure at top - pressure at bottom 
= (friction factor * downcomer length 
f downcomer diameter + entrance factor 

+ bend factor + exit factor) 

2 

* (downcomer mass flow rate) t (2 * (downcomer 

2 

cross sectional area) * density of fluid 
in downcomer * g c ) - density of fluid in 
downcomer * gravitational acceleration * 
height of downcomer t g 
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In the model notation 



PCPD - (fCD & Ldo/oCD +- EtOTftCO + 6&vJ0C0 

-pEXrreo) * M<LC *. * 2.0 /(2.0 *.*CJD*:X 2.0 
*. RB-CCJD ^ SC)— (eKCCO-^ S ^ SCO /<$C 



P6PK - (f 6H * LGtt/06tt EEnJT^SH + SEiCdSH 

+ EXfT6H')^ MSS * * 2.0/(2.0 ASH **2.0 
>fc RHOSH- * SC.) -KB06H * 6 * ^<Gh/gc 



for the screen and main bank downcomers respectively. FCD and 
FGH are friction factors of the form 



where R is the pipe radius and KS is the relative sand rough- 
ness 111] . 



F&H - 4*0/ (l.l^ >K AU3£ 10 (ksGH)) IiC;MeJ 



3 . Riser Pressure Drop 

The momentum equation for the riser boiling section 
must take two phase flow effects into account because the 
flow in the boiling section is not homogeneous. Thus, there is 
a relative velocity between the liquid and vapor phases here. 
The steam separators are included in the riser length. However, 



4- - i / (i.74 -z LjO$ 
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the effective length of the separators and thus the pressure 
drop are considered negligible because of a general lack of 
information concerning these items. 

The pressure drop due to single phase or homogeneous 
flow in the nonboiling riser section can be written 

SPF acceleration friction gravity 

^acceJora-f.'on ~ ^ 14 ) 

5 s * hrir jZ 

~ *4 % £ * lenflH) ^ ^ 

<&* *> * ISoJr^^ ^ 

2 / *. %<L 

Experiments conducted by Babcock and Wilcox [10] 
indicate that these homogeneous flow pressure drops may be 
modified to give the appropriate two phase flow pressure drops 
using correction factors that are functions of outlet quality 
and operating pressure. The two phase flow pressure drop can 
be written; 

AP = AP * r 

TPF acceleration acceleration 

+ AP * r 

friction friction 

+ AP . * r • . . 

gravity gravity 

where the r terms are two phase flow correction factors 
available by using curve fits of the data from reference [10] . 
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The form of these r terms become: 



f?S«WE - 2.4/144 * XFF^ ** 2 .^-b.SjzW % XFF0 + ,477fc l rc l 
kseavs. = 24.?14 * xu.^ * m z. f -i.sptic x xui^zl + .4 rile fxcj 

RJVa-E - ir.H^4xxFFjzi^ % Z.<p -H5.4444 *XFF0 ~.4<f40 r 7 [icj 
<AclT = 1^45-64^ X 142 S**. Z.f + IS.4144 *XU-0 -.$^7 |jo3 

KRICH = -34. * XFfjz5 a%z.(p +23. 7lb4 * XFF0 

4- . 5734 ^ 

i?Fi?iCK - -34. ^ x+qi 4: X 2.^ 423. 7lfe4 * XU^S frej 

+ .5734 



Therefore, in the boiling region of the riser 

4 p TPF= ^ r ad + ^-P 

S 6 ^ ^ (c\; 3 »VN©ier') 

4. <^. of ko>l.‘^ rgg.on)^,^ H^rdv/ 




The total pressure drop across the length of the riser is; 



AP AP spf + AP TPF 



In model notation the equations are 



38 



POPF = 



Prpi. 



(mcO * ;* Z.0 * (rhoec -RWODo)) / (< 5C % RHc££ * 
RWoOO * AOE * * 2.0) 3-(m.0 * FOE. * LOE 
34 MOO 34 * 2 . 0 ) / (S C 3K ODE * (rHOOD +• 
RHoee) AOE * * 2. 0) + (s *ed£ 34 (rhcod 
3-RHOEE.))/ (6C *2.0) +■ (m00 34 3< 2.0 34 
RAC.Le)/(sC * RHcEE.* AEF * *2.0) +- 
(2.0 * FEF * UEF* MOO * * 2.0 .34 RFRICE.)/ 
(SC * OEF * RHceE. * AEF * * 2.0) +- 
(<S * HEF * RHOEE * rsrav/e) / SC 



( pvnr * * 2.0 * (rhokx- rhott))/ <5c * rhck*. * 
Rl+OTT* A-JK* % 2.0)+- (“3.0* FTk^MST^ 
2.0 + LTk)/ (0C * DTK. 34 (RHCSTT4t?tJoKK)^ 
ATK*.* 2.0 + Q *£TIC*(RHoJT+-RHoKa)/ 

( 6 C * 2.0) +-((V1JT* 34 2.0 & RAClK.)/ 

(SC *. RI40ICK. * AkU * * 2.0) +- (2.0 34 
FKl * lkl* mts“**2.0 34 ^Fikick.)/ 
(sc * DK1_ * RftoKK. *. Aku * -34 2.0) +- 
(s * sku* AHokk * esieAv/K.)/Sc 
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The friction factors are in the same form as those 



for the downcomers 

foe- 1.56/(1/74 -z.4> * AU5SI0 (ksde)) trcl 

FEF - i. <pl (l. 74 - Z (/> * alosi^ ( <SEF) ) Crc 3 

Rnc = 1. 0 /(1.74 AloSI^OcSTK"}) \xc3 

FKu = 1.56/(1.74 -2.56 X Al£<Si 0 (K5ku)) \rcT) 



At a specified steady state operating point the 
pressure drop across the downcomers must equal the pressure 
drop across the risers and the downcomer flow rate must equal 
the riser flow rate. Therefore by equating the appropriate 
pressure drop equations and solving the resultant relation for the 
flow rates, the initial flow rates may be computed as 

KFF00 = ((kv$CD0 * G -*^C0 -G sK DO0 

^i^AV£)/((FC03K LC0/dCD -V-EnTOCD 4- 
GEiODCO 4 hxjtco)/(z.^^.aC£)^^ 2.^ 

^ -RHODO^')/ 

a Aof * )) +* 

(4.5^ ^ FOE. LDE^$ * 2 * 0 ) /(z .0 % DDF 

* (*HoE£^ +RHODO^) tK aof # 5K 2 , 0 ) +• 

RACL&/ (RH-o£E^ * AOF* * 2.^) 4 {4.0 
^ FEF LEF0 ^ &F£|CE.)/(z.^ ^ 00 F 

* R4oee^ * adf * # z.0))) % * <f>.s [rc] 



40 



N\u_0^= 

R&WWJK)/ ((F6H Ik 1.614 /O&H 4- E*Jr£6>H +• 

gekjosh 4- exitshVC*.#* A6H *: * z -0 

-V-^KCKK^ - RHo3TT$)/ 

(iRHOKK^ & KHCXJT0 *. ATL % Zffi)) +■ 

(4-5^ * FTK *. LJKjzi * z_fi)/(2.& X, OJ-L 

^ («HoliK0 +• RHOJTJzO & ATi- & * 2'<z0 +- 
RKH.Z-/ (RFOKKjzJ % ATL. * * Zj£)-t-(4.jzf 

*• FkL * LKL 0 * fcf=Rlck)/(z.jZS *orK 

* SHCKK0 * ACT. * * 2.0 ))) * * ^.s- ^ 



4 . Riser Continuity 

The relationship between the riser inlet and outlet 
mass flow rates is written in terms of the continuity 
equation for one dimensional unsteady flow, 



T e\l»\ - 




In model notation, this becomes 

IAll - rv\TT- DRBOTU % VOUTL- 
lr^FF =• fADD- DRPoOP & VOUDP 



l>e] 



The numerical differentiation technique used by CSMP- 
III is highly suspect, as are other techniques. An "averaging" 
system is used in the actual dynamic model. In addition, the 
flow rate down the downcomer and into the riser is held con- 
stant for the open loop boiler model. 
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5 . Riser Quality and Density 

The average density in the risers must be solved for 
explicitly. Linearly varying quality along the riser tube 
length follows directly from the assumption of uniform heat 
flux along the riser tube length. The average density along 
the tube length is the sum of the total change in density and 
the density entering divided by the riser tube length I 1 ] . 
Assuming the density varies only in the boiling length of the 
riser, the total change in density can be written; 



Since pC&) varies linearly from entering to exiting, p(&) can 
be written in the form 



The average density is 





where p^ = density of saturated water 
L = total tube length 





L-& 



where y = specific volume entering 
^out = < l ualit y at r i- ser outlet 
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Lg = boiling length 



LNg = nonboiling length 

V = change in specific volume from saturated 
° liquid to saturated vapor 

The integral of the explicit equation for p (£) may 

be solved after rearrangement to 

JL. 

(\/^ - XqJ: Lt^g v/j.^ + XoJ- 4 

Le ' L e 

which is of the form \ 

a + 

This yields the average density, 

p av * _L r — — i 

^ »- \X*Jt 

+■ bigg 

In model notation the formula is written 

KvcoFf ~ /((***#)* 

#. Alo£ (((XFFjz5)/VP)^VFei +1. 

•Vfcttooo^ * tottf) IxCjMej 

(LKi0/((XLL^)^VF€))#N 
A\_o^ (( XLL^ }K /VF +- 

>fc. lxk^S 



XqJ- x 1 

\k 
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To permit calculation of the quality term the energy 



balance on the liquid in the riser tubes is evaluated. Thus, 

rate of change = flow energy - flow energy 
of riser energy in out 

+ thermal energy in 

fL (jlW) - ~ +- <j,V, 

dX 

If the enthalpy term is treated as the average enthalpy in 
the riser the equation may be written 

In model notation this is written in the form: 

OHXLL- ((V\66 Q<q -Mll 

# XU, HFs/*.^) / (fcHOTL*. \ZDLKU} 

HKLL - \UT6,fU- ( h*LU^D\4XU_) 

- (hxlu - hf)* 

CMXFF - (MC£, (UOO-HF-XFF^ + 

-rv\FF * XFF tK FF$/z.^)/(fcl40DF *VOl£f) [f^Ej 
HXFF = UOT^euCHXFF^jOBXFF) [Vie] 

XFF - (v-VXFF- HF)^^.^/HRS pv»e] 

The initial conditions, HXLLO and HXFFO are calculated in the 
initial condition program. 

The steady state equations for quality are derived 
differently. The energy transfer rate to the riser liquid 
may be written in the form 

Q ~ A A, V\ 
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where 



AK — V\ oui* h.v> 

- - k; 0 

Solving for X yields 
out 

XoJ- ’ (q (kf - )) / (m 

or 

XFF0 - (6)1 4 f I^FF(ZS * (HF -H0 D^))/(m F F?S * HF6v) [xcl 
X Ll^ = (6)5# +• MU.0 * (rtF - HXJ0 )) / (FILLjZS * HF<S) [ r( Q 

6 . Riser Boiling Boundary Location 

The nonboiling length of riser tube is the product 
of the riser tube length and sensible to total heat ratio [1] 

1-WiS. — i. 

/ c it 

The sensible to total heat ratio may be expressed as the 
enthalpy change in raising the water to saturated conditions 
divided by the total enthalpy change, i.e. 

9 s - kj ~ 

4- ) — h.Vi 

The equations for nonboiling length are 

-LDF AmAXl( r (HF-VlDD^) ) 0.(^)/^(^oo^ 
t- XFFjZ $ ^ RF§) - HDO^ ) \ xej 

LTK56 - IXL^ AwAXI^RF 

i-XuL^ ^ HF$) - HvT3$) \Tc} 
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The Fortran function AMAX1 is used here because the 



above equations are used in an iterative loop in the initial 
condition program and there is a possibility during iterations 
of reaching a situation where the enthalpy entering is greater 
than saturation enthalpy. The final initial condition solu- 
tion prevents this. The AMAX1 function is not used for the 
dynamic model. Thus, 

L0£ = LDF as. (rF -Hdo)/((HF fXFF * HFS,)- HOo) [me] 
UTK = LTL * (HF - HXr)/(C HF +-*LL *. I4F$) -HTT) 

The boiling volume is derived directly from the solution for 
the nonboiling length. 

boiling volume = total volume * boiling length Hence, 

total length 

volef = voldf ^ lef/ldf t Me Q 

VOLK.L = VOLn. ^ LkL/_TL 

7 . Steam Drum Liquid Mass and Energy Balance 

The rate of change of liquid mass within the drum is 
equal to the sum of the mass flow rates entering and leaving 
the drum. The liquid mass flowing into the drum is the 
saturated liquid from the risers, the liquid from steam con- 
densation in the drum, and the incoming feedwater. The liquid 
leaving is that leaving via the downcomers. The dynamic model 
equation is 
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domsu - rau_ + mff * (i.^-xff) 

+ rv\CotJJ3 -+- W6jg — MCC - M<S €i 



The instantaneous mass of liquid in the drum is 

DmftSL ^ lOTSf^L. (orvi^Su^^PD^SL.) \j* e] 

where 

= (vbL.D'&V) ^ ^00&vn)/2,(2i [xc] 

The drum energy balance is derived similarly t i.e., 

ODWOHL — irALU * (1 .^-Xul)^ HP -S-^FFtK (l.^ - XFJ=) 
FF ■+• WCokJO H-F 6 -Fp^vSiS H 6 G — bACC, 
Yv\$$ H6.S ]poeJ 



The instantaneous drum energy is 



omortL - iiJr6^u(omoHL95 ) ornm+i.) 



I^eJ 



The initial drum energy is the product of the initial 
drum liquid mass and initial drum enthalpy. 



OnqOHL^ — DrvjASL^ ^ HD€u?vq^ XrcJJ 

The initial drum liquid enthalpy is 

DH ~ OrADHL ^ot*\f\ZL 

In the dynamic model, drum liquid enthalpy is considered the 
enthalpy of saturated liquid. 
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The steam condensation rate equation is based on the 
difference between the pressure and temperature of the steam 
and that of the liquid [12] 



M Co too - ^3 1 /(jmrA jzf ) ¥. % 

- tfSftT / (TS AT + Ht<zS. <j> ) * *. <f>. S') +- $ Z S& S' 




The water level in a marine boiler is generally con- 
sidered in reference to a mid-drum zero, that is, a water 
level of plus one inch implies the water level is one inch 
above the drum centerline. For computational efficiency the 
assumption is made that for small changes in water level 
around the midpoint of the drum the water surface area remains 
constant. This allows a simplified level equation, i.e., 

LE^/EU - (0rA0V- VOLOfcm ) / (LSTW\Oyv\ ^ 0 STp\ 0 m) JViej 

DMOV is the equation for the total volume of "liquid" in the 
drum. Recalling that a small percentage of steam leaving the 
separators is entrained in the liquid, the rate equation for 
"liquid" volume in the drum is 

DDrnoV - ((rv\FF + rv^LL-^,(3«s) % VF +Pdu. >K (mll. -Wvjff) 

W -V & VtSS — (fACCs ^ VF 
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where DMDVO is the volume occupied at time zero which is half 
the drum volume. 

8 . Circulation System Initial Condition Iteration 
Procedure 

The initial conditions for the circulation system 
are found by flow rate balancing in the downcomer-riser flow 
loops. An initial guess of flow rate is made using the assumed 
riser exit quality of .05 percent which is reasonable for this 
type of boiler [7], Coupling this assumption with the assump- 
tion that the percent carryunder is zero, an initial approxi- 
mation of the flow rate can be determined using the figures 
for initial heat transfer rate to the risers calculated in 
Section III-A. Hence, 



mass flow rate 



heat transfer rate 

assumed quality * latent heat of 

vaporization 



A first approximation of downcomer enthalpy for both 



banks can be calculated. This downcomer enthalpy is assumed 



to be the same for all downcomers. Using an energy balance 



on the liquid in the steam drum 

MFP^ - Q -l^/^XASurviE. 3^V4F6i) 

rv\LL^ - §5^/ C^A-Su^E ^ UPS) IXC] 

HCO^ = ((M,FF^ +rv\LL^- H-F -F * BeejzO / 

(MFF0-f^LL0) jxej 

- WCOtf LXCJ 



The screen riser inlet enthalpy is assumed to be equal to 
that of the screen downcomer; however the main bank fluid 
absorbs additional energy from the desuperheater located in 
that circuit. Therefore, the main bank riser inlet enthalpy 



4 ? 



must be calculated separately. Thus, 

HTT^ - 464^ 4- £rc] 

The downcomer density must be calculated for use in the pres- 
sure drop calculations. Hence, 

VC 0^ - ((fAFF^ VF +* yv\<3<3^ ¥ / 







tx<U 




— / Vco^S 


Drc] 


J^l4o6il+^ 




txc] 



The riser outlet quality is calculated using the 
initial quality formulation previously developed. This com- 
putation is followed by the calculation of the riser non- 
boiling and boiling lengths. 

The average density of the risers is calculated along 
with the two phase flow multiplication factors for use in the 
flow rate/pressure drop calculation. 

An updated flow rate is now computed and compared 
with the first approximation. If it is within a specified 
error criteria calculation stops. If not the previous 
approximation is updated and the calculations continued with 
the new approximation. 

Upon completion of the flow rate balancing the initial 
flow rates in the downcomer/riser loops are known along with 
the riser outlet quality-. 

The steady state drum specific volume and density 
are computed as 
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\/D£.mu^ - \fd 00 4- PCu % \)\/ 

RoofcmL - l.p /\/or^L0 



The initial steam mass in the drum is then calculated 

0STrv\<^ - VcuO&m R.Ho\/ 

9 . Superheater Pressure Drop 

All entrance, head, and exit losses are considered 
negligible in the superheater compared to the frictional 
pressure drop. 

AP = | L. 

V&Merace. 

By lumping the constants the pressure drop equation is 

Prv\rv\ — i°fsjfvj — KckjI. ^ ^ ^ ^ y/RBcWfO0 

The following definitions apply here: 

KoO 3 = (pmrvv^ + pMKJ^)/^l4o^n4^-VRHoN3|O0) ^<3 
K)4o/v\kJ(^ - (jRHowrv^ 4- 

KojOA. RHCM»O0)/ rAm ^ 



10 . Steam Valve Equation 

The flow through the steam valve is considered 
directly proportional to the outlet pressure and valve 
opening, i.e., 
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m - (L * P * (\feWe. Cpe^;^) 

rv\rv\i^v\ J_LL‘ - PWfU -5K KOS34 ^ VAu/€L ^Vv\e] 



The constant is computed in the initial condition program. 

koroM - Mwm / (vau/e$ ^ pnim^) £af[ 



E. EQUATIONS OF STATE 

The equations of state listed below are used in both the 
initial condition and dynamic programs. With the exception 
of the subcooled specific volume equation used for the 
feedwater entering the drum, they are reasonably accurate 
in the 300-1500 psi range. The subcooled specific volume 
equation is accurate in the £>Oo “ 8 'CO psc. range. 

PSAT * E*r(( AU3S 1X1 

TSAT ■* EXp((.2£lS*| AL06 (psAt) + 4.77123)) K3 

HSAT" - EXP (0.2fc4iT2, * ALO<5 (psAt) + 4«4&70?) W 

4-1.7, te-^4 

*PSAT* * 2.0-4.214E-0?** 3-0 IX 

RHOF -(oZ-%-0.0 i7?l^T3AT + 1.132 E.-0fT 
* TSAT * * - G.7?6E-^ * T-S/VT 

3.0 tX 

V66 ~ . ^Ib004?$ - .0000^2,0 1 4k ^T66 
-V 3. C?STHg. — 0? ^ Ttse 2.0 
- * TKS8 =+• ^ 3.0 + l-4jZte\E-13 

^ T£>8^ 4.0 - U14?E-\(* 5fcTT3e** 

+• ?.O2>4E“20*TB(5 *^6>.0 



52 



IV. RESULTS AND CONCLUSIONS 



A. GENERAL 

A listing of the initial condition program and the dynamic 
boiler model program is given in Appendices A and B respec- 
tively. The initial condition program output data must be 
properly formatted for input to the CSMP-III dynamic model. 

In addition, because the dynamic model utilizes the liquid in 
the steam drum as a saturation state point from which all 
other state points are derived, two of the initial conditions 
must be modified to eliminate a discontinuity between the 
steady state program and the dynamic model. The equations 
involved and an explanation are given below. 



(voLo*m L rc 3 

= DMA£Ljz£ 14 

D/vLDtf-L ^ liJrGfcL^DHL^OOMDH-O 

DPOftSL =. 0*<] 

0^ - Orv^DBu/OfviA^L \>' £ Q 



As stated previously, the enthalpy of the liquid in the 
drum is considered saturation enthalpy. The dynamic model 
must begin with the rate of condensation in the drum (MCOND) 
as close to zero as possible. This is a natural steady state 
posit ion - drum liquid and drum steam both at the same pressure 
and temperature. In order to insure this to be the case 
DMASLO and/or DMDHLO must be modified such that the initial 
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drum enthalpy (DH) very closely approximates the enthalpy of 
saturated liquid corresponding to the drum steam temperature 
and pressure (PMM and TMM). This is easily facilitated by 
the use of the CALL DEBUG statement in the dynamic program 
With this procedure: 

1. The initial condition program is executed and the 
initial conditions formatted for CSMP-III use. 

2 . The dynamic model is executed for only a short run 
time, i.e. 5 seconds. 

3. Observing the DEBUG output from the model, DMASLO 
and/or DMDHLO are modified such that PSAT and TSAT 
equal PMM and TMM. 

4. The model is reexecuted for a short run to check, 
MCOND should be very small. 

5. Since DH is implicitly related to PMM and TMM, the 
procedure may have to be repeated a few times. The 
objective is to force as many of the "DERIVATIVE" 
terms in the DEBUG output to comparatively small 
figures as possible. 

The model as used operates satisfactorily with DMASLO = 
8002.2 and DMDHLO = 3.96163E 06. This forces the initial 
MCOND term to .35949. 

Because of the CSMP function DERIV used in the program, 
the model is extremely sensitive to integration time step. 
Numerical differentiation is not a desirable function to 
perform in a dynamic model; however attempts to explicitly 
differentiate the equations concerned failed. One solution 



54 



to the problem is to "smooth" the derivative function by 
averaging it over several timesteps. The model performed 
satisfactorily with a fixed-step integration procedure 
(Runge-Kutta) , an integration time step of .04 seconds, 
and averaging the derivative over sixteen time steps. 

B. OPEN LOOP RESPONSE 

The open loop response to a ten percent increase in 
throttle opening is shown in Figures 2-7. As expected, the 
response of the main bank circulation loop and screen bank 
circulation loop is different. This is in keeping with the 
different purposes of those two loops, steam generation and 
furnace screening respectively. The effect of the valve 
opening increase is barely noticeable in the screen circula- 
tion loop. 

The swell effect is not noticeable. Further conversa- 
tions with Mr. Paul Weitzel at Babcock and Wilcox indicate 
that the percentage of "carry under" is not a constant as it 

was treated in this program. One to two percent is a good 

* 

starting estimate for steady state; however during a transient 
"carry under" mass flow rate is computed by subtracting the amount 

of steam leaving the boiler from the amount of steam produced. 

This should always be a positive number. The vapor that does 
not leave the boiler is "carried under.” The following program 
will implement this "carry under": 
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Modify 

the section of the dynamic 
model titled "COMPUTE THE 
DRUM SPECIFIC VOLUME" to 
read : 

PROCEDURE CRYUND=F I LTR8 ( MMM ) 
IF(TIME . GT . 0 . 0 )G0 TO 45 
CRYUND=PCU* (MLL+MFF ) 

GO TO 46 



ENDPROCEDURE 

DDMDV= ( ( MFF+MLL-MBB ) *VF+CRYUND*VV . . . 

+MBB0 * VBB- ( MCC+MGG ) * VF+MCOND . . . 
*VFG-RISE*VV) 

Follow the above statements with the unmodified 
equations for DMDV. 



The previous procedure was not implemented in the present 
dynamic model. The quality formulation in the present model 
is apparently too simplified and this causes a shortfall in 
quality at the outlet of the screen riser bank. As a result of 
this shortfall, there are instances during transient operation 
when the model is not "producing" as much steam as it is "using". 
Mr. Weitzel suggests using twenty node finite difference 
approximation for quality. This could be done using an equation 
for quality such as: 
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CRYUND= ( MLL*XLL+MFF*XFF ) -MMM 

RISE1=DELAY ( 250 , RISTIM, CRYUND) 

RISE2=CRYUND 

IF ( VALUE. GT. . 51)GO TO 55 

RISE=RISE2 

GO TO 57 

RISE=RISE1 

CONTINUE 
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where 



X = quality 

Z = distance up the riser (ft) 
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= heat input to riser (BTU/s-ft) 

o 

A = cross sectional area of riser (ft ) 

3 

= density of saturated liquid (lbm/ft ) 

V in = velocity of liquid entering (ft/s) 
h. = latent heat of vaporization (BTU/lbm) 

Simultaneously an implicit equation for continuity could be 
applied to produce a balanced mass flow rate without assuming 
that the downcomer mass flow rate remains constant . 

C. CONCLUSIONS 

The model presented is not a finished model. Further 
research is needed to successfully implement the shrink and 
swell theories presented. When the previously mentioned 
difficulty with riser outlet quality is solved, the model 
should approximate the dynamics of a wide range of D-type 
marine boilers depending on the initial conditions supplied 
to the initial condition program. 

The initial condition program develops the initial condi- 
tions necessary for the detailed boiler model with the 
relatively scant data available from the boiler technical 
manual shown in Appendix C and a small amount of data from 
common engineering handbooks. Because of the relatively 
small amount of data required the initial condition program 
and/or the model can serve as a basis for a boiler condition 
monitoring system. 
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D. SUGGESTIONS FOR FURTHER RESEARCH 



It is evident that the model suffers by using the CSMP- 
III function DERIV. In order to facilitate the explicit 
solution of the continuity equation for the riser banks some 
form of differentiation is required, either numerical or 
explicit. An explicit solution is much more desirable from 
a stability standpoint. A more detailed investigation with 
fewer assumptions might produce the correct form for the 
explicit differentiation equation. It should be noted, 
however, that when using explicit differentiation the model 
must average the derivative over at least two time steps to 
avoid the creation of an algebraic loop. This is easily 
done with the same PROCEDURE format used to average the 
derivative in the current model . 

The use of small perturbation techniques to linearize the 
model and thus allow a state-space representation for multi- 
variable control development and analysis should be undertaken. 
There is a computer code available locally to facilitate this 
for CSMP models; however it does not accept the CSMP DERIV 
function and/or the averaging procedure. A means of bypassing 
this problem would result in a general D-type boiler model - 
applicable to a very wide range of boilers currently in use - 
which could be linearized about specific operating points. 

These individual models could then be used for the development 
and testing of multivariable control systems. 

Using locally available optimizing computer codes, an 
optimal D-type boiler design could be attempted with regards 
to boiler geometry. 
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A more detailed analytical and experimental investigation 
of current D-type marine boilers should be undertaken locally. 
Little research has been done in this area and that which has 
been done has often been based on either an incorrect physical 
model of a marine boiler or on a nuclear steam generating 
plant. A starting point could be the data analysis of common 
horizontal and vertical steam separators followed by an opti- 
mum design for these elements. The improper and/or mainten- 
ance of drum internals apparently grossly affects boiler 
operation and water level stability during rapid transients. 
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APPENDIX 



A 



C8/3C/7S 2C.55.il 



ILt: CONSTANT FORT® AN M 



NAVAL POSTGRADLATE SCHCCL 



this progra* calcllates t H e constants and initial conditions 

FOR A US NAVY O-'IYPE BOILER GIVEN THE INPUT DATA OBTAINED F 5 C M 
ThE 6CILER T ECHNICAL 'Afl'jll., S T £ A'-' TABLES, AN Q STANDARD ENGINEER- 
ING HANJ8CUNS. the INPUT CATA ARE E\~ER£Q IN NAHE LIST FORMAT . 
NAMELIST DESIGNATIONS AND INPUT MEMONICS FOLLOW 



INCONl: 

OPPNT OPERATI JG PD l NT < PE RCENT ) * 

TOTSTM TTTAL STEAM FLO«<LB/HR) 

SFHSTM SUP ERF tATER STEAv FLCW (LB/HR) 

OSHSTH DFSUPERHC ATSD STEAM FLOW (LB/HR) 

P CRUM CRUM PRESSURt (PSIG) 

SHOT SUPERHEATER CUTLET TEMi>ERAT']=E (CEG-FI 
SHOP SUPER FEMER OUTLET PRESSURE (PSIG) 

OSHCT DESUPERHEATER lUTLET TEMPERATURE (DEG-r) 

D SHOP DESUPERHEATER OUTLET PRESSURE (PSIG) 

ECONIT ECONOMIZE 3 FEE 0 INLET TEMPERATURE (DEG-F) 
tCONCT ECCNC^IZEF FEED OUTLET T E m P E R AT : j R E (DEG-F) 

A I® IMP AIR INLET TEMPERATURE TO REG I STcRS (DEG-F) 

CJLTRF CIL T E “ R t c ATU RE AT EJ'NER INLET (DEG-F) 

AIRFLO AIRFLOW RA~E (LB/HR) 

0 ILFLC OIL FlO/, RATE (LB/HR) 

XCSAI3 EXCESS AIR (OcRCEN”) * 

DRAFT DRAFT LJSS (IN-H20) 

TC-A3SC FLUE GAS TEMPER ATURF LEAVING SCREEN (DEG-F) 

TGASSH FLLE GAS -EMPERATUPE LEAVING SUPERHEATER (CSG-F) 

tgasmb flue gas temperature leaving *ain samk <de g-fi 

TC-ASEC FLUE uAS T c'^P ER ATUR E LEAVING E r ON DM l Z ER ( 3EG-F ) 

TSCRN SCREEN TUBE WALL TEFPERA T URE (DEG-F) 

TSPHTR S'JPERFEATER TUBE WALL TEMPERATURE (DEG-F) 

HfiFVCL HtA T FcLEASc (.< B7M/ HR /C'J FT FURNACE VOLUME) 

FHAS FURNACE HEAT Ad SORB T 1 ON ( K6TU/HR/ SO FT 3ACIANT HEAT 
ABSORBING AREA) 

NOTE: "PERCENT" VALUES ARE TO BE LCGGED AS DECIMALS. 
i SCREEN 

D T UaSC AVERAGE INSIDE DIAMETER OF SC° EEN TUBES (IN) 

LA V SC AVERAGE LENGTH OF SCREEN TUBE (FT) 

NTUBSC NUMBER )F SCREEN TUB£S 

RHASSC RADIANT H E 4 t ABSORBING AREA OF FJPNACE SCREEN (SO FT) 
MASSSC TOTAL WEIGHT OF SCREEN TUBES (LB) 



1SPH7R 

AREASH TOTAL SUPERHEA" ER HEAT T R ANS F ER ARES (SO FT) 
MASSSh TOTAL WEIGHT OF SUPERHEATER TUbES 



SMNBANK 



OTUBM B AVERAGE INSIDE 01AMET=R OF MAIN BANK TUBES (IN) 
LA VMS AVERAGE LENGTH CF MAIN BANK TUBE (IN) 

NTU3M9 NUMBER OG -IAIN BANK TUBES 

MASS M B TOTAL WEICnT no MAIN BANK TU°ES (LB) 

AREAMB TOTAL HcAT TRANSFER AREA ( SQ FT) 



JECON 

DTLBEC .INSIDE DIAMETER CF ECONOM I 2E R T’J BES ( IN ) 
NPASSE NUMBER OF TUBE PASSES 
N T U SEC NUMBER OF t UBES °ER PASS 
LTUBEC AVERAGE LENGTH OF ECONOMIZER TIJ6E (FT) 
MASSEC TOTAL WcIGnT OF ECONOMIZER (LB) 



CONOOO 10 
CONOOO 20 
CONCCCAJ 
C3N00C40 
CONOOO 50 
CON CCCfcO 
CON 000 70 
CONOOO BO 
C3NCCGS0 
CONOOIOO 
CONOOI 10 
CON 00 1 20 
CONDO I 3D 
C0NCC140 
CONCG 150 
CGN00160 
C0NCC17C 
CONDO 130 
CONCOI 90 
CON CO 2 CO 
C0N002 10 
CONCC22C 
COM 00230 
CCNCC240 
CON CC2 50 
C3M00260 
C0N00270 
CON C02 30 
CO NO 0 2 GO 
CONCC3CO 
CON 00 2 10 
CQN00320 
C0NCC330 
CON JO j 40 
C0NC03 50 
CONOCO 60 
C0N00370 
C0NCC3 SO 
CON 00 2 90 
CQN004C0 
COMCCA 10 
C0NC0420 
CONGO A 30 
CON CCA AC 
C0N0045D 
C0NC0460 
COM CCA 70 
CONOOAdO 
C 0NC0A9Q 
CON 00 500 
C0NCC5 10 
CON CC 520 
C0N00530 
CONOC5AO 
CON C C 5 50 
CGNOOSoO 
C0NCC5 70 
CONDO 5 FO 
C3N00590 
C0MCC6CC 
CQN00610 
CDN00620 
CONCG630 
C0N006A0 
C0NCC65C 
CON C0660 
C0N00670 

con cceeu 

C0NCG6S0 

C0NC0700 
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oor>ponoononnooononononnonooor>nonnooo< >oo‘ 'nonoononnnnnoon^ >nnnn‘ n >no‘ >o 



FILE: 



CONSTANT FORTRAN Pi 



NAVAL POSTGRAOLATE SCHOOL 



$ C ES ° F 

OTUROS 

NTU8DS 

NPASSO 

LTUBCS 

APEACS 

MASSOS 



S ORMOCR 

DTUBDD 

LAVGC 

NTL6CC 



SHORDCR 

OTU3HC 
LA VHO 
NTU8H0 



SBCILER 

0 STM DM 

LSTMQM 

OWTRDM 

L WT R D M 

HMZRM 

HHOR 

hwTRCM 

FLRVCL 

NCTE: 

STHERMO 

H 5HQIT 
HDSOUT 
HECIN 
HECOUT 
< 1 )KF 20 

( 1 ) FFH 2 C 

( 1 )VSr.H20 
(2 )KSTM 

( 2) PR$TM 
(2IVSC3TM 

R 5HCU 
RDS3UT 
RFLUE 

NOTES: 



SLCSSES 
KSCSC 
KSDMB 
KS £CC 
K SOHO 
ENT SC 
EENCSC 
ENTMB 
BEN CM? 
EFTOG 
REND 00 

enthc 



INSIDE 3 1 AME TER OF DESUPERHEATER TUBE (IN) 

NUMBER IF OES'J°ERHE MER TU^ES PER PASS 
NUMBER CF DESUPERHEATER’ TUBE PASSES 
LENGTH IF 0E5U D ER HE A TF ~ TLSE (FT) 

T: T AL FEAT TRANSFER AREA OF QESUPERHE ATE* (SO FT) 
TOTAL HEIGHT CF DE S L PERM E AT ER ASSEMBLY (LB) 



AVERAGE DIAMETER CF DRUM OCWNCOME R TUBES (IN) 
AVERAGE LENGTH OF DRUM OOWNCQMER TUBE (FT) 
NUMBER OF CRUM CCWNCOMER TUBES 



AVERAGE INSIDE DIAMETER HEADER DOWNCOMER (IN) 
AVERAGE LENGTH CF HEADER C^WNCOMER (FT) 

NUMBER OF HEADER DONCOME^S 



DIAMETER OF STEAM DRUM (IN) 

L ENGT F OF S t EAM DRUM (FT ) 

DIAMETER CF WATER D FU M (IN) 

LENGTH JF WATER DRU W (FT) 

HEIGHT IF NORMAL WATER LEVEL ABOVE BENC.H mark (FT) 
HEIGHT CF HEADER (SCREEN) ABOVE BENCH mark (FT) 
FEIGhT OF *ATER ORU m ABOVE BENCH MARK (FT) 

FLFNACE VOLUME ( CU FT) 

CHOICE OF BENCH MARK. IS ARBITRARY 



ECONOMIZER CUTLET ENTHALPY (BTU/LBM) 

DESUPt RHEA TER OUTLET ENTHALPY (3TU/LBM) 

EN t HALPY CF ECONOMIZER FEED INLET ( B TU/ LB M ) 

ENTHALPY OF ECONOMIZER F c ED OUTLET (ETU/LEM) 

THERMAL CONDUC TIVITY JF WATER ( d TU /HR . F T-OEG-F ) 

P R AND T L NUMBER FCR WATER 

KINEMATIC VISCOSITY FOR WATER (LB^/FT-SEC) 

T H ER M A L CONDUCT IV IT Y FOR STEAM ( BTU/hR .FT .DEG-F ) 

PR AND 0 L NUMBER FOR STEAM 

KINEMATIC VISCOSITY FDR STEAM { LB M/FT— SEC ) 

SUPERhEA T ER 0U T LET DENS I T Y (LE'VCJ FT) 

□ESUPERHE ATER OUTLET DENSITY (LBM/CU FT) 

CtNSXTY OF FLUE GAS AT TGASSC (LBM/CU FT) 

(1) EVALUATED AT AVERAGE ECONOMIZER WATER TE M PE R A T U C E 

(2) EVALUATED AT AVERAGE DESUPERHEATER STEAM TEMPERATE 

ROUGHNESS RATION (SAND EQUIVALENT) FOR SCREEN TUBES 
ROUGHNESS RATIO (SAND EQUIVALENT) OF MAIN BANK ^UPE 
ROUGHNESS RATIO (SAND EQUIVALENT) F SR D R L M DDWNCO m Ep 
ROUGHNESS RATIO (SaND EQUIVALENT) FCR Hi EA CER DOWNCOMER 
SCREEN TUBE ENTRANCE LOSS 
5 CREEN TUBE BE^O LCSS 
MAIN BANK ENTRANCE LOSS 
MAIN BANK BEND LOSS FACTOR 
DRUM CCWNCCMcR E N T R ANCE LCSS FACTOR 
DRUM DOWNCOMER BE ND LOSS FACTOR 
HEADER CCWNCOMER ENTRANCE L~>SS FACTOR 



C0NC07IU 
CONGO 7 20 
CGNQ07 30 
C0NCC7A0 
CONGO 7 5G 
CON C C 760 
CUNC077C 
C0N0G7S0 
C0NCC19C 
CO N 0 0 8 0 0 
CON COS 10 
CON CG 820 
CO NOG 8 30 
CONGO 6 AO 
CON GO 8 50 
C0MCC860 
C0NCC870 
C0N008 30 
C0NCD8SU 
CON 00V 00 
CCN009 10 
CO N 009 20 
CON 0 C 9 30 
CGNQ0940 
C0NCC950 
CCNG0960 
CONG 09 7 3 
CON CG 9 80 
C0N0099Q 
CONG 1000 
CON CIO 10 
CQN01020 
CO N C 1 C 20 
CON 0 1 0 AO 
CONG 10 50 
CON C IC60 
CO NO 10 70 
COMOL 080 
CDNC1G90 
C0N01 100 
CONG 1 1 10 
CONCH 20 
CGN01130 
C0NC1140 
CON C 1 1 50 
CO NO 1 160 
CON C i 170 
CONO 1180 
CONG 1 190 
CONO 1200 
C0NC1210 
CONO 1 2 ZU 
CONO 12 30 
CON C 1 2 AO 
CONG 12 50 
CONO 12 60 
CON C 12 10 
CONO 12 30 
CQNG1290 
CONC 13C0 
CONO 13 10 
C0NC1320 
CONG 13 30 
CONCUR 
CON C 13 50 
CON 01360 
CON01370 
CON 0 13 SO 
CONO 1 390 
CONC 1 AGO 
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FILE: CONST ANT FORT 5 AN Pi 



NAVAL PGSTGRADLATE SCHCCL 



c 

c 



c 

c 

c 



C 



SENOriO HEADER DCUNCCMER 9ENO L^SS FACTOR 

IMPLICI’ 3 E A l ( A-Z ) 

INTEGER FPASS5 

DIMENSION ThZJI ZG)»Q°ASS(2CJ 
DEFINE NAMELISTS 



FAMELIST/ l NCC NO /CFPN7 , T OT STMjSPHSTM.CSHST" »PiDRUM,ShCT,SHOP»OSF r T 
IDSFOP. ECDJ l TtECON )F, MRTM? ,-JILFLO, AI RFLC, XCSAI R,DR AFT ,TGASSC . 
2TCASSF,~ CAS v» t r c-ASEC , TSC PM, TSPHT 3 , H 3 FVOL, FHA S, OIL IMP 
NAME LI ST /SCREEN/ CT’J3 SC , LAV SC , NTU 8S C » PH AS SC, M AS SS C 
NAM EL 1ST /’iPHT 3 /APf\SHt •’ A S S SH 

FA ME LIST/ v,N t JNX/ LT JB v 9 f LAV M3 , »j tij-jm a f M ASSM 3, ARtAMB 

NA ME LI ST /EC C'l/JTL DEC ,N PASS E . NTUfcEC , L TI Ofc C , MAS S EC 

NAM EL 1ST / D ES 3 F/ DTU BD 3 » NT’JR DS » "IP - SSO . L TU90S , ARc AD S » w A S SO S 

FA ME LI ST /J h m DC R / LT’J P CD » LAV CO * MU a DD 

NAMEL IST/hDROCR /DTUSHO ,LA VHD ..NTUBHD 

NAMELIST/ SOIL ER/ DST MCM, LST ''J'O OwTR DM, L WTR DM , HMOR M , HHQR , HWTRO M , 
SFURVCL 

NAMEL IST/THERMT/HS CTUT.HCSDUT, H.SCIN .HECDUT, 

1 KH2C ,PRH2C,'/; 02 C , <S T m, p r<; T‘i , v S FS T m , N SH'iMT . ROSO'JT , RFLUE 
NAMEL I sT/nSSES/ixSDSC ,KSO MB , isSJOD , KS CrlJ .ENTS C , EENCS C , ENT MB, 

$ 8 E N C M 6 » EOT CD. BE I CCJ, FN T HU . BEND HD 
FAMELiST / INCCiOl / ” RRO ,TVV3 ,TWWO , 

CTYYC.TZZC, IXXJ.Ti'D 
NAMEL IST/ l N C C N2 / "• A 1 0 . m F OOO » m AO 00 » 



CMNNO.MMMC.MQOO . 

CMC DO »M GHO » 10. 

GMFFO ,MILJ ,MM£, 

GMASSV.MA SS* .MASSY , 

CM ASS X* TASSZ, massor, 

CDMUHIC.JMASLD, JST v D ,CSTMOM, 

CL STM DM , DST I DM 

NAMELIS T / l-NCONS/ f vo, CAB.CPV, 

CCPY ,CPX ,CPZ , 

CCPW.COR .CRS ,CST . 

CCTU .CFF.G.GC 

FA ME LI ST/C00ST1/KSJL .KSDF . KSGH , 

CKSCC.KSCE.KSEF, 

CKSjK ,KSKl ,KZ , 

SKW.KPS.KSf, 

CXQM . KO'l-1 . KCN4 , 

C K X . KH J . K \ . 

CKTU, KY 

FAMELIST / CONS T2/ AJL, CJL.LJL2, 

CLJL1.DJF ,LDF , 

CLCFI.ACF.LCF, 

CAGH .CGil.DCO . 

CLC D , LDE 0 , L J< 0 

FA*EL IS" /CONST 3/ ACC,ZRENC1,Z3ENC2, 

C VOLJL » V C L ) F .VCLOhM, 

CVOLHJ.ZOF.Z JL , 

CLJL ,ODF ( OKL , 

GCEF.DJK.ZCD 

FAMELIST /CUNST 4/ F CO, FGH, FDE, 

GFEF ,F JK ,F XL , 

*ENTRGH,E ITRCD.BEODGH, 

C 3EFCCD, EX rjH, EX IT CC 

NAME LI S f /C OOST3/FHV , XLLO.XFFO , 

GTAM0, C50CFG, DRDJLU, 

GSIGMAA 

NAMELIST /OUTPL I / OIO.C •AO, 06C.06C, 09 0, MFFO.MLLO, CCO , 
GXFFO .XL LO, 3 FCUFO , CHOJLJ, RH'OFFO , 3 FiOLLO, 

GHSA T »HCDC »HAAO »H EBO .FiJJO ,HFG , RFOF, OM DVD , HX FFO , hXLLO 



/ Pc AO <5 , IACCNO ) 
READ ( 5, SCREEN) 
PEtCibtS'hT* ) 

* REA015 , MV3ANK) 



CONOID 10 
C ONO 1 420 
C0NC143C 
C0NG1440 
CONO 14 50 
CON C146C 
CONOID 70 
C0N0143O 
CJNC 1453 
C ONO 1500 
CGN01510 
C ON 0 1 5 20 
CCN01500 
CON C 1540 
CONG 15 50 
CONO 1560 
CON C157G 
CO NO 1580 
C 0 N C 1590 
CO.NClcCO 
CONG 16 10 
CON C 1 620 
C3N01630 
C0NC1640 
CON C 16 50 
C0N01660 
C0N0167C 
CON 016 80 
CC NO 1 6 90 
C0NC17CG 
CONC 1710 
CO NO 17 20 
CON C173G 
CONO 1 740 
CONC 17 50 
CON C 1 7 tC 
CONO 17 70 
CON Cl 7 EG 
CON 0 1790 
CON Cl 8 00 
CONC 18 10 
CONO 18 20 
C0N0183C 
CONC 18 40 
CQNC1850 
CONC 1860 
CON C 18 70 
CONC 1 8 80 
CON C 1 8 50 
C0N01900 
CONO 19 10 
CONC 1920 
CONO 19 30 
CON 0 19 40 
C0N0195O 
CONC 1960 
CON C 1 9 70 
CONO 1 980 
CON C 1 5 50 
G3N0 2000 
CON C2 010 
CON C 20 20 
CO NO 2 030 
C0N02C4C 
C0NC2050 
CCN02060 
CONC 20 70 
C0MC2080 
CONO 2090 
CON C2 ICO 
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FILE: CONST£N T FORTRAN PI 



NAVAL POSTGRADUATE SCHOOL 



C 

c 



\ 



RE ADI 5.ECCN) 

REAC (5. C 65° t- J 
RE AO (5 » 0 F J JCi) 

REACH 5, H OR DC 3 ) 

PE AO (5 , e n 1L f> ) 

READ ( 5 • THE nV 0) • 

REAO( 5*LOSjcS ) 

CALCUIAT1CA CF P R EL I N IN APY CONSTANTS 

Pl=3.1A15927 

G=32.2 

PR AA=PRH21 

V I$CCA=VSCF23 
THCONA=KH2 J 
PR AN- PR STN 

VI SCCN = VSCS T ^ 

THC CNN = K E T N* 

XASU" E=0 .0 5 
FC L=. 01 5 

GC = 2 2 .2 

VALVEO=OPP'JT 

KSCO=KSDHD 

ks OF=K S DSC 

KSDE = KSOF 

KSEF =KSDF 

CM=. 11 

CPV=CP 

CP t=CM 

CPX= C M 

CPZ=CN 

CPW=CN 

K3GH = KS OCO 

KS JL=K SD F9 

KS JK=KSJL 

KS K l = KS J L 

ENTRCO=E NTHD 

ENTFGH= ENTOC 

9ENQCD=9EI\QHD 

R ENOCH* '3ENJM-1 

EX ITCD = 1 .0 

E XI TGH = i. 0 

PS AT = PCRU M + 1A .7 

HSAT=EXP{0.26A52*aLCC< PSA’ ) +A .46 703) 
HFG= S22 . 15- C. AC5lc -P SA T+ 1. 717F-OA* PS 
J-V=FSAT -h-FG 
HMMO*H V 
HF=HSAT 

TSAT = EXP ( .22151*ALCu (°S AT J+A.7 7123 ) 

MF CQ C=0 1 LF LO /3<->G C. J 

A* A C C 0= A HFL-I/ 26uk..O 

TRBO= TjA S SC 

TVVO=TSCRN 

TVW0=TSPFTR 

TQRO = TRR C 

TSSO=TGASSH 

THKO=TS AT 

TTTO = TGAS-*e 

TliUO=T GASEC 

TTLO= ( TTTO + TUUOJ /2 .0 

TA M< 3= 8C . C 

“MMO=T IT STm /3600 .0 

MMf'E=M"MC*.0D01 

K'MMIII=SPHST"/j600.0 

TNNO =SnC T 

TP P C=D SHOT 

PNN0=SF1 3 +1 A. 7 

Ff»KO=FORLH 

MMM I I=JSH3TN/3oOC.O 

MNNO=Hyvi I 

TAAQsEC CM T 



C0N021 10 
C0N02120 
C3NC2 130 
C0NG2L AO 
CONC2150 
C0NC2U0 
C0N021 70 
C3NC21 60 
CON 02 l 90 
C0NG22 00 
COM C 22 10 
CQN02220 
C0N022J0 
CON 0 2 2 AO 
CQN022 50 
CON C 2260 
CONC2270 
CO NO 2280 
CON C 22SO 
C CN02300 
COMO 23 10 
C0NC2320 
CGN023 30 
CON C23AC 
COM C23 50 
C0MC2360 
C0.NC23 70 
CQN02360 
C0M0239O 
CON 0 2 AGO 
CONO2A10 
CON C2 A 20 
C0NC2A30 
C0N02A AO 
CON C2A5C 
C0N02A6G 
CON02A70 
C0NC2AS0 
C0N02A90 
CON C2 5 CO 
CO M 0 2 5 10 
C0MC2520 
CON C 25 30 

■'** 2 . 0-A. 219E-09*PSAT**3 .0 CCK025AO 

CO NO 2 5 50 
CON 0 2 560 
C0N02570 
CON C2 560 
C0NG2590 
CGN02600 
CON C 26 10 
C0N02620 
CONC2620 
CON 026 AO 
C0N02650 
CON C266G 
C0NG2570 
C0N026S0 
CON C26SO 
CONJ2700 
C0N027 10 
CON0272O 
C0N0273O 
CON C2 7 AO 
C0N027 50 
C0NQ2760 
CON 027 70 
CONO 2 780 
CONC279C 
C0NC2SCQ 
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Fi IE : CCf^STANT FCP r R AN 



naval POSTGRADUATE SCHOOL 



n 



TR BO-EC 0 ^OT 

V860= .01600488- .0 000 0 20 146 *T8 R C+. 00000003651 1* TRBO**2.0 
£-8. 142c- 1 L* TBoO** 3 . 0 + l . 408 1 E- 1 3* T 6 00 * *4 .0-1 .14 3 E- 16 * 

6TBR0 **5 .'J>o.u34E- i J*T*B u** 6. C 
CT CTO=n PF7 CL* IQG 0 .O/JoOO .0 
HNNQ=H5HCUT 

C FU R 1 6= F H as * 1000.0/ 3C00.0 
A R A C=RH A S $ C 
HP P C- HI) SOL T 
L C F1 = L AV SC 
L0F2=LDEl*NTUBSC 
CDF=0TU3$C*NTU3$C/12. 0 
CQ E = CCf 
OE F=OUE 

ACF=f>l*CTU3SC**2.C*NTUBSC/<4.C*144. 0) 

VC lDF*ACF*LuFl 
' MASSV=M A $3 SC 
M A $S V = M A SSS H 
LJL1 =LA V.*3 
LJL2*NT»JBM3*l JL 1 
o ji=DTi.ja^a*N T udM e/ iz .o 

D JK=D JL 
CKL=DJK 

AJt- FI* CT* I .O^NTUB^B/ (4 .0 <144.0 )• 

VJLJL-A JL^L JL 1 

M ASSY= M ^SS‘1 3 

►A J S*=MASSfcC 

M A S SZ = i A S 50 S 

Z B ENC1 = H lORN-HHOo 

LDF = Z EE NCI 

ZDF=Z Q ENC 1 

Z 6 EM C2 = h CR f*-H^T POM 

L J L= ZBE NC2 

ZJL=ZBENC2 

ZCO=ZCF 

ZGH = Z J L 

CROC F0 = 0 .0 

ORCJL 0=0.0 

LAR=LTUREC 

CAe= ctusec/ u.o 

AR EAEC = fsTU3EC*LA6*t 3 I -GAB* NPAS S E 

ON P= CTURCS / 12.0 

IC0= LAVHC *N T U 0HO 

LG H = LA *00*NTUBDD 

DCC*OTUBHO^ITUSHC/ 12 .0 

DGH=LTuecc-‘jruaoc/i2 .0 

ACD=NTU3H J* ( M* ( C TUB HO /2 4. 0 ) **2.0 ) 

AGh=MTUeCC*(° l*( C r U c 00/ 24. u ) **2.0> 

VOLD FM=( PI *03 TMOM** 2.0/ (4.0*144 .0) ) *LST M C M 
VOLHJ=( p I *O wTR’Jt4**2. 0/< 4. 0*144) ) *LWTROM 
HA AO =HEC IN 
H3 8 0=HEC CLT 

RHQF = 6 3 .tf-U. 0173 1*T$AT+ 1. 132E- C5*TSA T**2. 0-6. 7S6E- 03* TS AT**3 . 0 
VF =1 .0 / PhCF 
VFG*S24. C/PSA T-0. 1 
. VV=VF*VFG 
RHO V =1 .O/VV 
RHOFG=RHOF-P H9 V 
PHQM *0 = R h CV 
RHONNO=R SH3 L T 
RH0P°0=R CS9UT 
M A AO = MN*0 
M3 3 C=M A A C 

C CALCLLATE THE TCTAL MASS CF FLUE GAS IN THE FURNACE 

C 

. MA$SCR=RFLUE*PJP.VOL 

C N 

C CALCULATE THE TOTAL ENERGY ENTERING BOILER 



CCN02810 
C ON C 28 20 
C0NC2E30 
C0N028 40 
CON C2E50 
CQN02860 
C0NC2370 
CO N C 2 d 80 
CQN02890 
C0NQ2900 
CON 0 29 10 
CON02920 
CON C 29 30 
C0NC2940 
C0N02950 
CON C2960 
CGN02970 
CON C 29 80 
C0NC2S90 
CQN03000 
C0NC3C10 
CON 0 30 20 
CQNC3030 
C0NC3C40 
CON 030 50 
CO NO 30 60 
CON 030 70 
C3N030 30 
CQNC3090 
COM 0 3 1 CO 
CO NO 3 1 10 
CONC312G 
CON 03 130 
CON 03 1 40 
C0NC31 50 
CON C 3 1 60 
CGN031 70 
CON 0 3130 
C0N03190 
C0N03200 
CONC32 10 
CO NO 3220 
C ON C 32 30 
CON 03 240 
C0NC325J 
COM 03 2 60 
CONO 32 70 
C 0 MO 32 80 
COM 03 290 
C0N03300 
C0MC3310 
C0MC3320 
CCN03330 
C0NC334C 
CON 033 50 
C0NC33 60 
COM 033 70 
CQN03380 
C0NC3390 
CON 0 34 00 
CONC3410 
CON C 2 420 
C0N03430 
C OKO 3440 
CONC345U 
CQN034o0 
C0NC34 70 
CON 0 34 80 
C0N03490 
COM C3 5 CO 
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FILE: CONSTANT FORTRAN p I 



NAVAL POSTGRADUATE SCHCOL 



OQC=CTOTC a FURVCL 

CALCULATE " LUMP EC" LOWER FEAT IN G VALUE OF FUEL/AIR M I X 
F HV = C C 0/ M F 0 C 0 • 

CALCULATE ENERGY ABSORBED BY RADI AT I CN IN FURNACE-C 
Cl C=CFLF/e* ARAD 
C2 0=C10 

CALCULATE THE S T E F A N- B C L T Z F A N CONSTANT MULTIPLIED BY TFE 
SCREEN AREA 

S IGM AA=Q 10/C ( TQR C+46C. 0) ** 4. 0- (TVV 0+460. 0**4-. 0) 

CALCULATE THE TCTAL FLUE' GAS ‘MSS FLOW RATE INTO BOILER 

FCC0 = FFCC0 + *ACC0 

MRRU=MQQC 

f$S 0= m RRC 

FT T0 = F ESC 

MUU C=H TIC 

COMPUTE THE ENERGY TRANSFERRED TO THE WATER ILDWER) DRUM 
BY THE OfcSUPERHEATER 

Q90=FNNO*M NNO-H PPO ) 

099Q=C9O 

CALCULATE THE ENERGY TRANSFERRED TO t H E MAIN BANK RISERS 

C6C=FFF0* ( H. M M 0- Ft 6 E 0 ) -C2 0-099 0 
0 5 C= Q 6 G 

COMPUTE THE FAIN BANK HEAT TRANSFER CO EFF I C I ENT -TU BE M ET A L 
TO STEAM SIDES 

KY=AREAF3/i. 782EC6 

CCFFUTE THE M A I N BANK TUBE METAL TEMPERATURE 

TYYO= (060/ <KY*PSAT**< 4.C/3 .0 ) ) >**( 1.0/3. 0I + T3AT 

COMPUTE THE MAIN BANK HEAT TRANSFER COEFF I Cl ENT- FLUE GAS 
TO TUBE *ET£L 

T ST0=( TSSO+TTTO ) /2. C 

KST = C5 0/ C-'SS0**0 .6* (T STO-TYYO ) ) 

COMPUTE THE SCREEN BANK HEAT TRANSFER COEFFICIENT-TUBE 
FETAL TO STEAM SICES 

KV = C 1 0 / (PSAT**(4.0/3 .0 )*( TVVO-TSA T) **3. 0) 

COMPUTE THE HEAT TRANSFER TO THE SUPERHEATER 

04 C= M MMO* ( HNN3— H FFJ ) 

030=040 

COMPUTE THE SPECIFIC HEAT CF THE STEAM 
IN TFE SUP E 3 FEATER 

CMN=Q 40 /( MMMO*( TNNC-TM«0) ) 

COMPUTE THE SUPERHEATER HEAT TRANSFER C3E FF I C I ENT-FL J E 
GAS TO TUBE ' META L-C 

TRSC=(TRRC+TSSO) /2. C 



CON 035 10 
C3N03520 
COM C3 5 3C 
CON 03 540 
C0NC355O 
CONC3560 
CON 03 5 70 
CON035 80 
CONC35SC 
CON 03 600 
C0NC3610 
CON C 36 20 
CO N036 30 
CQN03640 
CON03650 
C0N03660 
C0NC2670 
C 0*1 C 3 6 SO 
COMO 3690 
CON C37CC 
CON 037 10 
CDNC3720 
CON C 37 30 
C0N03 740 
CONG 3 7 50 
CON C 3 760 
C0N03770 
CON C 3 7 SO 
C0NC37S0 
C0N03800 
CON C 36 l 0 
CON 0 3 320 
C0NCJ830 
CON C 3640 
C0N03350 
C0N03660 
C0N036 70 
CONO38B0 
CON C 28 SO 
C0NC39CO 
CO N C 3 9 1 0 
CON C3S2C 
CUN03930 
CON C3 9*+0 
CON C 39 50 
CCN03960 
C0N039 70 
COM 0 29 80 
CO NO 3 9 90 
C0NC4C00 
C3NC40 10 
CO NO 40 20 
C0NC4C30 
CON 04040 
C0NC4 050 
CON 040 oO 
C0N04070 
CONC4C80 
C0N04090 
CQN04L00 
CON C4 1 10 
C0NC412O 
CON 04 130 
C0NC4 1 4 C 
CO NO 4 1 50 
C0NC416C 
C0NG4 170 
C0NO4180 
C0NC4L9C 
CON 04200 
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FI IE : CONSTANT FORTRAN Pi 



NAVAL POSTGRADUATE SCHOOL 



KR S= 030/ ( M R R 0 * * 0 • 6*( TR SO-T *W0] ) 

COMPUTE THE HEAT TRANSFER TO THE ECONOMIZER 

C 8 C= N A A 0* ( H 6 60-h A A 0 ) 

Q 7 C= 0 8 0 



COMPLTE THE SPECIFIC HEA^ OF THE FEEDWATER 
IN THE ECONOMIZER 

CAB = Q8C /( *AAC* { T3 60-TAA0) ) 



CONFUTE THE SPECIFIC 


HEAT 


OF 


FLUE 


GAS 


I 


N 


THE 


FURNACE 


CQR= CCQO-QLO )/ (M PRO* (TOR Q- 


TA MB 


n 














COMPUTE THE SPECIFIC 


HEAT 


□ F 


FLUE 


GAS 


I 


N 


THE 


SUPERHEATER 


CRS=Q3G/(MRRC*(TRRC-TSS0) ) 


















COMPUTE THE SPEC IF I C 


HEAT 


OF 


FLUE 


GAS 


I 


N 


THE 


VAIN BANK 


CST=G5 0/ (MSSO*( TSSU-TTTO ) ) 






m 












COMPUTE THE SPECIFIC 


Ht AT 


OF 


FLUE 


GAS 


I 


N 


THE 


ECONOMIZER 



CTl=070/{ VTTQMTTTO-TUIJO) ) 

ccvfute the economize? heat transfer coefficient-tube 

METAL TO LI CUlu 

KX = ( £• 02 3^ TrCONA )/CAE J*(4 ,0/ < P I*0AB*V IS C0A*NT J 9EC U**C .8 
£*PRAA**C.4*A~ EAEC 

COMPLTE Trtt: L.MTC FOR T HE ECQNOM I 2E R 

LMTD A 6 = 08 3/ { KX *M 1 AO * *0 .8 ) 

COMPUTE THE ECONOMIZER TUBE METAL TEMPERATURE 

EXP0EC = E*P( { T69G-UAC) /LMTD4B) 

TXX0=<TAA0-TeRu^EX : >0 bC)/< I .u-EXPOGC ) 

COMMUTE THE ECONOMIZER HEAT TRANSFER COfcFFI Cl ENT-FLUE 
GAS TO ~UBE MgTAL 

KTU=Q70/{ M TTJ *( TTUO-TxXO) ) 

COMPUTE THE DESUPERHEATER HEAT TRANSFER CC£FF I C I EMT-T'J 9 E 
metal TO STEAM 



K Z=( ( .G23*TH03NN)/3NP>*< 4. 0/ C PI *0NP* V I SC QN*NTU3DS ) 3 
£**O.8*0F AN**0 .i* ARE ACS 

COMPUTE THE DESUPERHEATER LM TO 

LMTONP=QSO/(KZ*MNNC**C*3) 

COMPUTE THE SPECIFIC HEAT OF STEAM IN THE 0ESU° ERHEATER 

CNP=C90/ (MNN3*(TNNJ-TPpO ) ) 

COMPUTE DESUPERHEATER TUBE METAL TEMPERATURE 

E XPOD $=E * P { < TNNG-1PPC) /LMTCNP) 

T Z 20= (TNNJ- EXPO OS *T D P U ) / { 1 .O-E >P33 S) 

COMPUTE THE 0E SLPERHE ATER HEAT TRANSFER CGEFFICIENT- 
WATER DRUM LIOUIO TO TUBE METAL 



CO NO 42 10 
CON C422C 
CON 042 JO 
CON 0^2 AO 
COM C 4250 
CO NO 4 2 60 
C JN04270 
CON 042 80 
CQN04290 
CONC43CO 
CON 043 LO 
CON04320 
CONC433C 
CGN04340 
CON C4 35 0 
CONC4360 
CCN04370 
CON C 43 £ 0 
C0NC4390 
CJN04400 
COM 04 4 10 
C0N04420 
C0NC443G 
C ON04440 
CGN04450 
C0NC4460 
C0NC4470 
CQMC4480 
C3NC44SC 
CON 045 CO 
C0MC45 LO 
CON C45 20 
CON 045 30 
C0NC4540 
CON 045 50 
C0N04560 
C0NC45 70 
C0NC45E0 
C0NC < t590 
CON C46CC 
CON 046 LO 
CONC4620 
CON C 46 30 
CO NO 4 6 40 
C0NC465C 
CON 04660 
CON04670 
COM C4660 
C0MC46S0 
C0N04700 
CON 04 7 10 
CO NO 4 7 20 
CO NO 4 730 
CON C4 740 
C0N047 50 
C0NC476C 
C0N04770 
CON04780 
C0NC4790 
C0N04800 
C0NC4810 
C0N04320 
C0N0483 0 
CON 048 40 
CQNC485G 
C0NC486 0 
C0NC467C 
CON 048 60 
C0NC4890 
C0NC49C0 
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FILE: CONSTANT FCRTRAN PL 



NAVAL POSTGRADUATE SCHOOL 



KH J= C99C /( TZZC-TSAT) 

COMFIT E The SUPERHEATER L* T 0 

LMTCM^ (TANO-TMMO ) / ( *LDG ( (TWHO-T^mq j /( TWWO-TNN C) ) ) 

CALCULATE THE SUPERHEATER HEAT TRANSFER COEFFICIENT-TUBE 
M fcT £L TO ST £ AM 



KW=Q40/ ( *MMj**o. 8*L M TDMN ) 

CALCULATE THE THROTTLE VALVE FLCW CCEFF i C I ENT 



KGN4 


-MMM 


1 1 i /( 


VA 


LVE 0* FfsNO 


) 










CAL 


CUL AT 


E 


THE SUPEP 


HEAT ER 


n UT LET 


DENSITY 


K0N3 


s(3M 


MO+PNNO 


) /( 1HGMM0 


♦RHDNN 


C) 






PH CM 


(SO - ( 


RnC WM 


3 ♦ 


S hCM'lO )/ i 


.0 








KOM 


= i ( P 


FFL-F 


MOJ^HONO 


) / MM MO 




.0 






CCM 


fUTE 


TH 


c\ F R I CT n 


N F A C T 


CRS 


FOR 


THE 00 *N COM EPS 




ANO 


THE 


RI 


SE R S 










FC 0- 


1.0/ 


( 1 . 74 


-2 


. 0*ALC310 


KSCD) 


} 




m 


FGH= 


1.0/ 


( I. 74 


-2 


. C*\l GGIOK SGH ) 


) 






FOE = 


1.0/ 


(1.74 


-2 


.o*:lcgio 


( < s C E ) 


) 






FEF = 


1.0/ 


( 1. 74 


_ 2 


. C*AL JG10 


( KSfcF) 


) 






F JK= 


1.0/ 


(1.74 


-2 


.u*iL T GIG 


( i\S J K ) 


) 






FKL- 


1.0/ 


(1 .74 


-2 


.C*-W_CG13 


USkL) 


i 







START ITERATION TO BALANCE CIRCULATION LOOPS 



COMPUTE INITIAL VALUE OF RISER OUTLET MASS FLCW RATE 

MFFC=0 10 /( XA SUME *HFG ) 

ML LQ= C5 0 / (X ASUM E*FrG) 

CALCULATE THE INITIAL DOWNCOMER ENTHALPY 

7L HC0C = ( ( MFF0+MLL0-MB30 ) *HF+ MB eO *H8BO +0 .0 * (MFF3 + ML‘L0 ) 

£*HV )/( MFFO+MLLO) 

HGHO = H C CO 

COMPUTE THE INITIAL DRUM ENTHALPY 
HORUMO=HCOO 

COMPUTE MAIN BANK RISER INLE T ENTHALPY 

Mh FO - *LLO 

H J JC-HGH C + 09 0/MHH0 

CALCULATE THE RISER INLET DENSITY 

RHCOOO^RHOF 
* RH C J JO- PH CF 

COMPUTE THE COwNCOM ER DENSITY 

VCCC= ( ( MFF C+M LL O-M3 0 C ) * VF+ MB B0*V8 BO) 

£/ (MFFO^mllO ) 

RHOCOG^i . O/VCDO 
R HO G HO = R HOCOO 

\ CALCULATE THE RISER OUTLET QUALITY 



CON 049 10 
CON 049 20 
C3NC4920 
C3N04940 
CQNC4S5C 
CON 04960 
CCNC4970 
CON C 49 80 
CON C4 9 90 
C0N05000 
C3NC5C1C 
CON C 50 20 
CON05O3O 
CONC5C40 
C0N05050 
CON0506 0 
CON 05070 
C0N05080 
C0NC5C80 
COM C 51 CO 
COK051 LO 
CON C 5 i 2 C 
CO N05130 
CON C 5140 
CON C 5 l 50 
CO NO 5 1 60 
C0NC5I 7C 
C0N051S0 
CON 05 I 90 
CON C 5 2 CO 
C0N052 10 
C ON C 52 2 C 
CON C 52 20 
CONG 5240 
CON C 52 50 
CON 0 52 60 
C0N05270 
CON C5 2 80 
CDN05290 
C0NC53 00 
CON C 5 2 10 
C0N05320 
C 0 NO 53 30 
CONC524J 
CO NOS 3 50 
CON C 5 2 60 
CON 053 70 
CO NO 53 80 
CON 05 290 
CON05400 
C ON 054 10 
CON 0 54 20 
CO NO 5 4 30 
C0NC5440 
CON 054 SO 
C ONO 5460 
C0NC5470 
C ON C 54 60 
CQNQ5490 
C0NC55C0 
CO N 0 5 5 1 0 
C0NC5520 
CON C 55 20 
CQN055 40 
C ONO 55 50 
CON 055 60 
CO NO 5 5 70 
CON C 5 5 60 
CON C5 5 SO 
CON 05 6 00 
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FILE 



\ 



: CONSTANT FORTRAN P l 



NAVAL POST OR ACU AT E SCHOOL 



HDCO=HCCO 

XFF G-( Jl Cw'FFO* ( HF-HCOO ) ) / (’*FFC*HFG ) 

XL L0= (Q50*MLLo*( Hr-HJJ 0) ) / ( MLL C HFG ) 

CALCULATE THE RISER OUTLET DENSITY 
’ PHCFFC=SFCF-XFFO*PHPFG 

rhollc=phof-xllo*rhofg 

CALCULATE THE NCN30ILING LENGTH OF THE RISERS 

ZB EN C 1= FNORM-HHCP 
LOF = ZBE NCI 

L0E0 = L0F *4MAXl({HF-HCDC),0.0)/l (HODO+XFF 0*HF G) -HOD 0 ) 

ZBfcNC2=HNC P M — hWT £ CM 

LJL=ZPENC2 

LJK0=LJL*4MAXL( (HF-HJJO) ,0.0) / ((HJJO+XLL 0*HF G ) — H J J 0) 

CALCULATE BOILING LENGTH CF RISERS 

LEF0=LDF -LUE0 
LKL0=L JL-LJKO 

CALCULATE THE AVERAGE GENS IT Y IN*HE SCREEN RISER 

PH0CF0= (1 .0/LCF ) * (LEFO / ( (X FFO) *VFG )* ALOGI ( ( XFF C 
G) /\F)*VF3*1.0)*RHCUC0*lDEO) 

CALCULATE THE AVERAGE DENSITY IN THE MAIN BANK RISER 

RHCJLO* ( 1.0/LJL ) * ( L N L C / ( ( X L L 0 ) * VF G ) *AL0G( ( XLLO 
G * VFG ) / VF *■! « 0 ) *KriC..w0*LJK0 ) 

CALCULATE T HE EFFECTIVE HEIGHT OF THE RISERS 

ZDE0=Z0F-LEF0 
ZJKO = Z Jl-LKLO 
ZEF0=ZCF-Z0EC 
ZKL0=ZJL-ZJ.<0 

COMPUTE THE TWO PHASE F LOw MJ LT I FL I C ATI ON FACTORS 

PG PA VE=24. 794*XFf 3**2 .0-6.5066 *XFF0+ .97 76 
RGR4VK= 24. 794*XL LC**2. 0-6. 5 J6 6* XULO+ . 9 7 76 
RACLE=L5.4’564.*XFF0**2.G*15.49‘,4*XFFC-.OOCO7 
RACl J=l 5.4564«X LLO** 2.0+ 13 .4944*XLL0-. 90007 
RFR ICE = - J4.0S22-*XF=0**2.0 + 2 3. 7lo4* xFF 0+ .3734 
RFRICK=-34.03Z2*XLL0**2.0 + 23 . 7 I64*XL L0 + . 8 734 

CALCULATE S ECONO APPROXIMATION CF MASS FLO* RATE AT 

EX IT OH SCREEN P IS ERS M 

PHO E EO= R HO 000 

HF F CC= I (PHCCOO*G*2:C-G*ZCEO*( ( R HOD 00 ♦ RhOF EO ) /2 . 0 ) 

G-G*Z EF0*R FOES 0*R GRAVE ) / ( ( rCD* LCD/ JC J + EN“ RCG 
G+eEMCG + EX IT CO )/ ( 2.0 <490**,: 0) + ( I 0 H 16EO-RHODDC ) 

G/(RHOEEO*RHCODO*AOf** 2.0) ) + <4 .0*FGE*LDE9 *2 .0 
C )/ ( 2.0*D0F*I RH3£fcC+ c HODDO) *ADF**2. 0) +RACLE/ 

G <RHCEEO*Al'F**2 .0 ) + <4 .0*FEF*L EF0« FC ICE)/ 

G( 2. 0*DOF*RHQEEO**uF**2.O) ) )**0.3 

CALCULATE SECONC AP FR CX I^at I ON CF ”ASS FLOW RATE 

AT EXIT OF .MAIN BANK RISERS 

RHOKKC=PHO J JO 

MLLOO=( ( RhuGHO*G*ZGH-G*ZJ:<0*( ( RHOJ J 3+°H0KK0) /2. 0) 
G-G*ZKL3 *PHL’KKJ*A. v; PAV:< )/ l < F CF*L CF/OGH + E'I 7 = GH 
G+8ENQGH+E <1 IGH 1/(2. 0*4 Grt* -2 . 0* RH CGHO ) + ( ( FHUXKO -RHCJJO ) 

G/ (PHCK.K J * 3 h"lj JO *CJL **2.0 ) ) +( 4 . C*F JK *LJK 0*2. 0 
G ) /I2.C*0JL* (RnJKKO+RHQJJO ) *AJL**2.0 ) +RACLJ/ 



CON C 56 10 
CON 05 1 20 
COMO 5 6 30 
COM C 56 4 C 
COM C5650 
C ON 05 6 6 0 
COM C 56 70 
C0N056 30 
COM0569Q 
CON C 5 700 
CO NO 5 710 
C0NC57 20 
CON05730 
C0N0574O 
CONC575C 
COM 05 760 
C0NC57 70 
C0MC57E0 
COM 0 5790 
COM05800 
COM 0 58 10 
C0N0582O 
C3NC53 20 
COM C 5 6 40 
C0NQ58 50 
C0NC566C 
COM 05 8 70 
CONC58SO 
CON C5E90 
CON 05 9 00 
C0MC59 10 
C0NC59 20 
CCN05930 
C0NC5940 
COM C59 50 
C0M05960 
C0MC5S 70 
COM 0 59 80 
CCNC5990 
CON C60CO 
C0N060 10 
C 0N0602G 
COM 060 30 
C0N06040 
C0N06C50 
C0MCOC60 
CCN060 73 
COM C6 CSC 



CO M 06090 
CQNC61 00 
C0MC61 10 
CON C6 I 20 
C0N06I 30 
COM C 6 14 C 
C3N06I50 
C0NC6 160 
CO MC 6 I 70 
C0M06I33 
C0NC6ISC 
COM 06 2 00 
C0NC62 13 
CON C6 2 20 
C0N06230 
C0N06240 
COM 06 2 50 
CO NO 62 60 
COM C6 2 7C 
C0M0623J 
C0NC&290 



CON C63C0 
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FILE: CCNSTAM FORTRAN =1 



NAVAL POSTGRADUATE SGHCOL 



£1 RH1KK C*AJL**2. 0) + ( A. OF KL* LKL0*RF RICK) / 

£ ( 2.0 *CJK *Ph‘)K>i\ 0 *AJl_ *»2 .0 ) ) ) ** 0. 5 

COM PAR F PREVIOUS APPROXIMATION FCR RISER MASS 

FLO* C ATE TC CURRENT, IF WITHIN E°R0P CRITERIA CONTINUE 

IF NOT, UPDATE AND T hEN REITERATE— 

CHEC K = 0 .0 

IF IABSI MFFu-MFFCCJ .LI. .01 )G0 TC 52 
MFF0= (MFF0J--'FF0)/2.C+MFF0 
CH E C K = L « C 

62 IF(Ab5|MLLC-^LL0CJ.L7..O1)GO TC 54 
NLL0= l ML L )0-MLL0 ML .0+ M LL0 
GO TO 71 

64 IFICFECK .EQ .l.OJGC TC 71 

CQMPLTE INITIAL MASS OF LIQUID. IN CRUM 
DMA S 10 = ( ^CLDPM*ahCCCO )/2 .0 

COMPUTE ThE INITIAL DRUM "LIQUID" VOLUME 
0M0VC=V0LDRM/2.0 

CCMPLTE ThE INITIAL ENERGY STORED IN DRU 1 * LIQUIC 

CMChlO=CNASLO*hCPUMO 

EQUATE INITIAL FLOW RATES 

MCDO=MFFC 
MGHQ= M LL 0 

CQMPLTE INITIAL MASS CF STEAM IN STEAM 
CRUM 

OSTMO=V3LORM*RHOV/2. C 

CQHPLTE HXFFO AND HXLLO 



C 



6C0 



6 Ci 



\ 



hXFFO = HF 4 FFC-*XFFO /2 .0 
HXLLC=hF 4 -HFG*XLLC/ 2 . C 



WRITEIS 
FORMAT! 
WRIT E ( o 
WRIT E 16 
FO PM Ml 
WRIT E (6 
WR I T E ( 6 
WR ITElo 
WRIT E (6 
WRITE (6 
WR ITElo 
WRITE 16 
WRITE 16 
WR IT E (6 
WRI TE ( 6 
WR ITElo 
WRITE (o 
WR I TE l o 
WRIT Elo 
WRI TE 16 
WR ITEl 6 
WRITE 17 
WR I TE ( 7 
WRITE l 7 



,600 ) 
IHL) 

, IN CON I) 
,601 ) 

1 H 0 ) 

, IUCON2J 
,601 i 
, IN CON 3) 
, 601 ) 

,CQ 4ST I) 
, 60 1 ) 
,CCNST2) 
, 6CL ) 

, C0N3M ) 
,601 ) 
,C r JNST4) 
,60 I ) 

,C CN 3T5) 
, 60 l ) 
, 600 ) 
VlU'PUT) 
, IN CON 1 ) 

, I NC CN2 ) 
, INC ON 2) 



WR I T E 1 7 , C C NS T l ) 
WRI TE l 7 ,C ON ST 2 ) 
WRITE 17 , CON ST j ) 



CON 06 3 10 
CON 06 3 20 
CON 06 3 30 
C 0 N 063 AJ 
C 0 NC 6250 
CON 06360 
CONG 6 3 70 
CUNC 63 EQ 
CON 06 390 
C 0 NC 64 00 
C 0 NC 64 10 
C 0 N 0 o 420 
C 0 NQ 6430 
C 0 NGo 440 
CONO 64 50 
CON C 6 A 6 0 
C 0 NC 6470 
CQNC 6480 
CON C 64 SG 
CON Oo 5 GO 
C 0 N 065 IG 
CON 06 5 20 
CGN 06530 
C 0 N 065 A 0 
CGN 06550 
C 0 N 06560 
C 3 NC 6570 
C 0 NC 6590 
CON 06 590 
C 0 NC 66 CC 
CO N 0 6 6 1 0 
CON 06620 
CONC 6630 
C 0 N 066 A 0 
C 0 N 0665 Q 
CON 06 660 
CONO 6 6 70 
C 3 N 066 E 0 
C 0 NG 6690 
CONOo 7 00 
CONC 67 LG 
CO N 0 6 7 2 0 
CQNGo 730 
CON G 6 7 AO 
C 0 NC 6750 
CON 06 7 60 
CON C 6 770 
CO N 0 6 7 8 0 
C 0 NC 6790 
CONG 66 CO 
C 0 N 068 10 
CON C 6 E 2 C 
CO N C 6 0 3 C 
CQNC604O 
CON C 6 E 50 
C 0 N 06860 
CONO 60 70 

coNGoaao 

CO NO 6 090 
CON C 6 5 CO 
C 0 NC 69 L 0 
C 0 NG 6920 
CON C 6 S 3 G 
C 0 N 069 40 
C 0 NG 6950 
C 0 NC 6960 
C 0 N 06970 
CQNC 630 C 
CON C 6 990 
C 0 NC 7000 
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FILE 



CONSTANT FORTRAN P 1 



NAVAL PQSTORACLATE SCHCCl 



WRITE (7 f CCN" T 4i 
WRI TE ( 7 ?CONST5) 
WRITE (7, OUTPUT) 
STCP 
END 



CQNC7C1C 
CON C 7020 
CQNC7030 
CON C 70 AO 
C0N07050 
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APPENDIX B 



Cfc/2C/j« 21.26.1? 



FI IE: CSHF F CRT RAM Pi 



NAVAL PO 57GR S J LATE SCHOOL 



//HOPEFULL JD6 i 2473, C567.NE 31 ,,9) WALKED R y C 1319* ,TIY£=2 
II EXEC CSf / °XV 
//>. SYSIN 00 * 

* 



* 



* 



* THIS 0 7V°E RULER. MODEL IS WRITTEN IN CSYP-llI 

* C C y P ) T E 5 LANGJACc.- IT |S A N^ML INEAR M ID EL 

* WHICH UTILIZES A SET OF INITIAL CONDITIONS A NC 

* CT HEP DATA PREPARED BY AN INI T I A L CCNDIT1DN FORTRAN 

* MCCEl WUCh PRECEDES THIS ONE. ^ id £ M JDcL IS GOOC 

* FOR S y ALL TRANSIENTS AF DO NO TnE CPcRATiNG POINT 

* specified 3Y THE INITIAL cono it idn program. 



<$****4********** = S>**=:<'*'** = *iS:t *« ** :£« =* * * ** = *c #»**£***$:*:**$=**** = *:=$* 

INITIAL 

r ONS 7 DR A4 1 = C. C ,DRA42 = C. 3 ,G PA43 = 0. 0 , CRA^AxO.D , 0PA45=J .0, OR A46=0 .0 , . . . 
OR /4 7=0.0, C? A48= u .0, DR 121= C.O, OR,- 32= C. C, OR A 3 3=0. 0, DR A 34= C. C,. .. 

DRA35 = L. C ,L'RA36=J.0,DR;37=O.0 ,CR A33=0.0 
CONST ON E= 1. C, ... 

T RP,0= 2563 .0000 ,TVVJ= bOO.OuOCO , TOWO= 1033. OOOC ,TYYC=. 

515.42749 ,TZZO= 600. 99976 ,TXXO= 360.36450 
CCNST TWO* 2.0 * ... 

DHJHO = 2 754667.0 ,CHASSO= 5415 .3164 , y AAO= 45.00C030 ,TAAO=. 

2 54. C C C Go , y FOOO= 3.307',999 ,YACCO = 55. 396330 , MNNO= . . . 

7.7500000 , NMvtj* 45 .000000 ,VALVEQ= 0.50999999 ,MQQO=... 

55.2C6375 ,CVN= 0. {,5523023 ,CA3 = 1 .0761400 , C°V = 0 .10999995, .. 

C°Y= 0.10999995 ,r2x= C.1C599995 ,CPZ= 0. 10999995 ,C°W = 

0. 10 9999Q5 , D**JnLO=39ol62o.O, C’' AS LO =600 2 .2000 ,DMDVO= 163.65696 
CONST TH0£E=2.C, ... 

DRO CFO= 0.0 , D c ° JL 0= 0.0 ,M'1-1E = 0. 4499 9979E-02 , . . 

MGhD = 062.70325 ',MR a O = 45.000000 ,lST'in y = 

16.665 95 6 ,DS1MDM= 6G.COOOOO , y C 00 = 960. 02 734 , Cl CVO=l 63 . 66 

CCNST FO'i R = 4 .0, ... 

KSJL = C. 55999991 E- 04, hSDF = 0 .2999 9996 E-04 , KS Gh= 0 .4 9999 D99 E-05 , KSCD= . 
C. 49999999 E- 05, XSJE- C. 25999996E -04,KSEF= 0. 299999 96E-04 , ... 

KSJK=0.59999D91E-D4, ... 

KSK L= C. 5 5 59 9 5 9 IE -04 , A Jl= 10.673426 ,OJL= 190.57166 ,LJL2=... 

26420.000 , L JL 1= lC.COOOJO ,DDF = 2 5. 39 7 4 7 ci , i_OF= 1 5. 75 OOOC , . . 

IDF 1= 12.7CCC30 ,ACF= *.0646043 ,LGH= 25.00o000 ,AGh=... 

1.4451332, SAN 3 L=250,RI STI M* 10. 0 
CCNST F I V E=5 .0, ... 

DGH= 1.9183321 ,DCD= 4.6699991 ,ICC= 9-,. 500000 ,LDEO = 

2.6166315 , L JK 0= C. bC 5C0 102 ,ASD = 2 . 3547 77 3 , -‘A S S V= 7951 . 3 934 . 

,MAS 5 W = 10752.000 , V ASSY= 26192 .000 ,Hl5 J x= 2444 2. 0 JO, M A S3Z =. 

1617. CC C C , Z3ENC 1= 1 5.750C00 ,t-CO= 0 .3 10 23 69 3E-0 1 , ... 

FCc= 0 .927143465-01, FEF= C .927 1 4646E- Oi ,F GH = C. 8 1023 6 93E- C 1 
CCNST SIX =6.0, ... 

Z BENC 2= 12.41CCC1 ,3ASSDR= 71. 355 ,VDLJL= 1 06 . 73426 , VC LC . 

F= 51.620463 , VTLDR v = 327.JU96 ,VOLHJ= 96.174744 

CONST SEVEN = 7.0, ... 

COR = 0.31599614 ,CFS= C. 3045367C ,Z ST= C.290423J4 ,CNJ = 

C. 26 026 064 ,CNP= 0.533471 39 ,KMP = 16.508d20 , K W= 2 . 0997 72 4 , . 

S IGMA A= C. 2 1 132812E-C9 ,RR S= 0. 87139436 ,KST= 2.2 037563 

CCNST E1GFT=6.0, ... 

KTU = C. 441266CC , EN TP GH = 0 . 99999 964 E -0 1 , E NT PCD= 0. 999 9 9 964E- 0 1 , . . . 

BENDC0= C . 99994964E -0 1, EX 17GH= 1 . 0 0000 CO , E < I TC D= . . . 

I .0000000 , hAAO= 221.00000 , K.DN 1= 0. 1064p8 7bF- C 1 , . . . 

KDN 4= . 1C 777515 ,F JN= 0.961 96 149E-01 , F<1 = 0 .9 6 1 Do 1 *9 E-0 1 , . . . 

GC= 32 . 1999 9 7 ,TAM8= 80. CCCJCu ,ZCD= 15. 75 0000 ,ZDF* 

15.750300 ,ZJL= 12.410001 ,LJL= 12.410001 ,... 

3 ENO GH= 0. l49»55Sd,NCNJ = t4 7. 71460 ,G=3 2 . 1 99 9 97 , 0 ST YD =249 .59802 
CCNST N1N£=9.0, ... 

KX= 1C. 03 73 73 ,KZ= 14 .930669 , KV= 0 . 32 769 o l 3E- 0 5 ,.. . 

FHV= 193CC. 496 ,0DE= 35.39747b ,DKL= 196. 5 7166 ,OEF=... 

35.697476 ,CJK= 19b . 5 7 1 06 ,N Y= 0. 266 6 2 50 2E- 02, . . . 

ORA 1 C=C. C ,0RA2 0 = 0.0 ,XF F D = 9 . 1395 JO 00 E -01 , . . . 



CS*'000 10 
CSYC0020 
CSMCCC20 
CS “00040 
CS y CCC5C 
CS -, OOOtO 
CS v 000 70 
C S YCCCEO 
CSMoC090 
CSV00100 
CSMCC110 
CS'-OO 120 
CS Y00130 
C SMOO 140 
CS Y00150 
CS M CC16C 
CS“00 170 
CS YOO 180 
CSMCC190 
CS y C 02 CO 
C S v 002 10 
C SHCC220 
CS '••00230 
. CSY00240 
CSMC0250 
CS YC02 60 
..C S M C C2 7C 
CS yOO 280 
CSYCC290 
. CS y CC3 CO 
. . CS yoo 3 10 
C SMCC22C 
CS M 00330 
. CSMC0340 
. C SM C C 3 50 
CS y 00360 
CSY00370 
..CSM00380 
CS YOO 3 90 
CSYCC400 
CS y 004 10 
. CSM00420 
C S V C0430 
CS v J0440 
CSYC0450 
..CSN00460 
.. CS YOO 4 70- 
..C SMCC4 60 
CS ,J, 00490 
C SYCC5C0 
C S y C 0 5 10 
. .CSY00520 
CS M CC530 
CS YU 0540 
..CS YC0550 
. . C SMCC560 
CSM00570 
CS Y00580 
CSMCC55C 
CS Y00600 
C SYC06 10 
CS ‘*006 20 
, . CS Y00630 
C SMCC64C 
CSYC0650 
C S V CG660 
CSMCC670 
CSM0J680 
CSYC0690 
CSM00700 
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f i l£ : C S V P FORTRAN PI 



NAVAL POSTGRADUATE SCHOOL 



XLLO= C.4CC (CL CCE-C l,MFFC=St C. 0273-*, Ml L 0=662. 70925 
CONST T^io .0 , . . . 

OR A1 1=C. C,DR A 12 = 0.0 ,0RA13 = C.C f DRA14=C.0,DPil5=0.0,CRAlo=0.0,... 
DR A1 7=0.0 i OP A 2 1= C. o * D : 4 22=0.0,0* A 23 = 0. 0, DR A2 4 = 0. 0 . . 

FCL=. 015»CC 0=63836. 3 95 , l GF = 12 .4 1 0001 , .. . 

QR A2 f = C. C, 0 = 4 26= C. 0, Uft A 2 f=C. 0 »riX LLC = 505. 31 ,HXFF0 = 502. 36 
CYNAH1C 

4 

* 1MPUT EQUATIONS 

* 

VA LVE=. 0154PAVPI1 3.0) — .015 *RAMP(20.0 1+.51 

m\!N=MNNC 

VAA=MAAO 

1AA= TA AO 

* 

* CCMFUTE The TOTAL FLUE GAS FL^W RATE INTO BOILER 

* 

KC=“CCO 

* 

* COMPUTE THE ENERGY ENTERING THE BOILER 



CO =Q00 



4 

* 

* 

* 

a 

* 

4 

4 



CO w F'JT E the ENERGY TRANSFERRED TO THE SCREEN RISERS 
VI A RADIATION 

=SIGVAA* ({TRP+460 .0 ) **4 .O-ITVV +460 .01*44.0) 

COHP'JTE ThE RATE EQUATION FOR FLRNACE FLUE GAS 
Tc V FePATURE 



D T PR=(CQ-C1-'1RR4C0R4(TPR-T A^3 1 )/( M 4SSQR*CQR) 
* 



* COMPUTE THE FLUE GAS TEMPERATURE 

4 



TRR=IN TGRL ( TRRO ,DTRR) 

* 

4 CC M F'J" E T F£ T E*‘FcR4TJRE OF THE FLUE GAS 

4 LEAVING THE SUP fcR HEATfcR 

4 

PHI 1=2. C4N'RP=4C.4 = C«S/KPS 

TSS= ( TPR*(PHI 1- 1. 0) +2. C*TUU) n PHI n-l .0) 

4 

4 CONFUTE ThE SUFEPhEATER ENERGY TPANSFER 

4 FLUE GAS TO T'JP 6 METAL 

4 



Q3=MRR*CRS*< TRR-TSS1 

4 

4 CC v FU‘ r E The '“AIN BANK ENERGY TRANSFER 

4 FLUE GAS TJ TUBE METAL 

4 



Q5=M«S4C ST4 ( TSS-TTT) 

4 

4 CO^R’ITE THE TEMPERATURE OF THE FLUE GAS 

4 LEAVING THfc VAI N BANK 

4 

CH 12 =2 .0 4'iSS * 4Q .4 4CST/ KST 
TTT=(TSS*{PHl2-1.0)+2. 0 *TY Y ) / ( PH 1 2 4-1 .01 

4 CCMFUTE TFE TEMPERATURE OF THE FLUE GAS LEAVING- 

4 THE ECONOMIZER 

4 



PH13=2.C4C TL/ML 

TUU= (TTT4(RHI2-l.C)+2. C « T X > J /(PHI 3* 1.0) 

4 C 3 M PL YE ECCNC M I iE P ENERGY TRANSFER 

4 FLUE GAS T) T'JPE METAL 

4 

Q7=MTT*C TL* < TTT-TUU) 

4) 



CS M 007 10 
CSMCC720 
CSM 00 7 30 
CS M C0740 
CSMCC750 
CS M 00760 
CS-00770 
C SM00 780 
CS M 007 90 
C SMC03 00 
CSM009I0 
CS. MOOS 20 
C S v C C630 
C SM 003 40 
CSPC0850 
C S‘-' CC 660 
CSM00370 

c s^cosao 

C SM 00390 
CS y 00900 
CSM0C510 
CSMCC920 
CSM009 30 
CSMCCS40 
CS M 009 50 
CS VC0960 
C SMCC570 
CSM0U980 
CSMQC990 
C SMC 1000 
CSMC1010 
C SMC1020 
C SM C 1C 30 
C S m 0 1 0 40 
C S^Cl C50 
C3«0 1060 
CG F C1070 
CS-'CICEO 
C S v C 1090 
CS NO l 100 
C SM c 11 1 0 
CS-'J 1 120 
C S MCI 130 
CSM C 1 140 
CS mo 1 1 50 
C S*C11 to 
CS k ‘0 1 170 
CS MCI 180 
C SMC 1 ISO 
CS" 01200 
C s*- 012 10 
CSMO 12 20 
CS ”01230 
C SMC1240 
C SM C 1 2 50 
CS V 01260 
CSMC1270 
CS M C1280 
CS V Cl 2 90 
C S" C12C0 
CSM013 10 
C S VO 1320 
C S M 0 13 20 
CS V01340 
C? M C1250 
CSMQ12C0 
CS 1 3 70 
C SM C 1 3 EC 
CSMO 1390 
C S W G1 4 00 
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FILE: CS* d 



FIR TRAN PI 



NAVAL POSTGRACUATE SCHOOL 



♦ 

* 

* 

A 



ISE CONTINUITY re lati cnshi f to compute THE MASS 
FLOW 3 ATE UP TH E 5 I SE P S 



L= M J J-DRHOJL*VOLJL 
PFf=' , CC-CPFOCF*V r, LCF 
ORA1 E = ut F l W J-Al 0 ,RnOJLJ 
U 3 A2 8=CE 3 I V ( OP A 2 0 t ? H 00 F ) 

PRCCECLSt CRFC JL = F IL T F 1 (DR A14) 

IFIKEEP.NE .1 ) G C TO i 

DPFDJL = ( C 3 A 10 + CR A 17 + DR A Lfc+J= A I 5 -pOR A I 4+D 3 A 13+DR A 12+DRa l 1+ . . . 

DRA31+LR AJ2 + CRA33 +CR A34+CP A35 +CRA36 +D c A37+D 3 A3d ) / 16.0 
DR A53=0R A52 
DRA52=CRA51 
0RA51 = C.R 431 
CR A31=CR 132 
CR 432= CRA33 
QRA33=DRA34 
CR A3 4=CR A3 5 
DR A3 5=CR A36 
DR A3t=DRA 37 
CR 43 7=CR438 
DR A3 S=DRA1 1 
CRA11= CR 112 
DRA12=CRA13 
DR A 1 3=DRA 1 A 
CR A 14= CP A 15 
DkAI 5=DRA16 
ORA 16=0R A 1 7 
DRA1 7= CR A1 8 

1 CCNTl.NLE 
tNCPR^CEDURF 

FFCCECURE CPFCCF=FIL~R2(DRA24) 

IFJKFEP.NE. 1)00 TO 2 

CRFOCF=( CR A2S + CRA27 0R A2o + C° A2 5+D 3 A2 4+.9RA 23+D R A 22*0R A2 1+ . . . 

0 R A4 1 + CR A42 + CR A43 +CR A4a + GR A45 +D.RA4u +D P A4 7 +DAA43 ) / 16 -0 
CR A 6 3= C 5 A 6 2 
CR 46 2= CR A6 1 
DRA61=DRA41 
CR A4 1= CR A42 
0RA42 = CR 443 
DR A43=0R A 44 
CR 44 <*= CR 445 
DK A45=DRA46 
C® A4 6= CR A 47 
DRA4 7=CRA4 8 
DR A 4 6= DR A 2 1 
CP A2 1=CP A22 
DRA22=0RA23 
DRA23=DR A24 
CRA24=CRA25 
DRA25=DR A26 
C° A 2 E= OR A 2 7 
CP 42 4= CR A2 8 

2 C ON TIN LE 
ENCPROCtCURE 
* 

* COMPUTE THE AVEPAGE DENS1 T Y IN TFE RISERS 



RHOJL=( 1.0/Ljl>*( (LKL/(XLL*VFG))*AL3G(( X L L / VF ) *V FG *1 .0 » . . . 

♦ PhC J J * L J K. I 

RHCDF = (l .0 /LCF J* ( (LEP/ (XFF*VFG ) ) *ALOG ( (X FF/VF )*VFG*l .0 » . . . 

♦ PH1CC*L J£ J 
RHCFF==FCF-XFF*PF1FG 
RHOL L-RHOF-XlL^RHOFG 

* TOMPUrt THE MAI i RANK ENcPGY TRANSFER - 

* TueE H E t A L TC *4IN BANK M 1XTU 9 E 



* 

C6=KY*FS 4T**U.O/2.C)*( TYY-TSAT) **3.0 



C SM C141C 
CS" 0 1420 
C S MCI 4 30 
C SM C 1440 
CS HQ 1 450 
c S M C 1460 
: 3" 01470 
C3 v 0 1480 
CSMC14S0 
CS M C 1 too 
CS"0 1510 
CSMC1520 
CS MO 1 530 
CSM01540 
3 SMC 15 50 
CSMC1560 
C SM01570 
CS“01530 
CS V C1590 
C $*'016 00 
C3 "0 16 10 
CS v 0 1620 
C SM C 1630 
CS MO 16 50 
CS"C 16 50 
Cd"Clt60 
CS MJ 16 70 
: $ " C 1 6 8 0 



C S*'Cl 7 00 

: smc 17 io 

CS M J 17 20 
CS MCI 730 
C3 M 01 740 
C 3MC1 750 
C SMC 1760 
C3MJ1770 
C S V G l 780 
C3"0 1790 
CS M 01 800 
CS M C1610 
CS *•0 1 8 20 
CS V C 1 8 30 
C S 1 C 1 0 40 
CS MO 1 8 50 
G SMC i860. 
C 3 ■* C 1 8 70 
CS •'0 18 80 

c s" c i Fsa 

CS" 01900 
C S" 01 9 10 
C SMC 19 20 
CS MO 19 30 
C S MO l 540 
CS-O 19 50 
CS MO 1960 
C S ' ' C 1 9 7 0 
C S" C 1 9 80 
CS *0 1990 
C S' C2LC0 
CS" 0 20 10 



C SMC2 020 
C SMC2C 10 
CSM02040 
C S M 0 2 0 50 
09*02060 
CS V J20 /O 
CSML2C80 
CS M 02090 

cs v o<.ioo 
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PILE: CSMP 



FORTRAN PI 



NAVAL POSTORACLATE SCHCCl 



* COMPUTE THE SCREEN RISER ENERGY TRANSFER - 

* TUEE METAL TO SCREEN RISER MIXTURE 

* 

02=K V*PS AT **( A.C/3.C)*! 7VV-7SAT) ** 3. C 

* COMPETE THE RISER AVERAGE ENtHALPYS 

* 

FDF=HOD*XFF*hFG/2 .0 

FJL=HJ J+XLL*HFG/2. 0 

* 

* CCMFLTE THE RISER OUTLET ENTHALPYS 

* 

hFF=FF*X FF*FFC- 
FLl=HF-t-XLL*HFG 

* 

* COMPUTE THE RISER EFFECTIVE HEIGHTS 

* 

IC E= ZCf-L E F 
Z Vi« = Z Jl-LKl 
ZKL=ZJL-ZJK 
ZEF=ZCF-ZDE 
€ 

* COMPUTE THE SPECIFIC VOLUME OF THE ECONOMIZER 

* OUTLET LIQUID 

* 

V R B= .0 160CAeo- -C0C0020 1A6*73B+. 00000003651 l* T» B**2 . 0. . . 

-8.1 A2E-11 * T 6 3**3 -) + i. i0it-l.J*T3B«*A .0- 1 . 146 E- 16*=T3B **5. 0. . . 

* e. 0 2AE- 2 C*TB8 **6. 0 

4 

* CCMFLTE The RATE tOIIATICN FOR RISER OUTLET QUALITY 

* 

ChX L L= (MGG*( FJ J-FF-XLL *FFG/ 2.0 ) +Qo- MLL*XL L *HFG /2. 0) / (RHO JL* VOLKL) 
HXLL=1 NTGRL IHXLLC .GhXLLl 
aLL= ( H >L L-hF I * * . 0 /HFG 

CFXFF=(MCC*( l-CC-FF-XFF*FFG/2 .0 )+Q2-" FF*X FF* FFG /2 .0 ) / ( R HOOF * VOLEF ) 
HXFF=! NTGRL ( 3XFF0,DH<FF) 

XFF=(FXFF-FFJX2.J/HFG 

«r 

* 

* 

* CCMPLTE ThE RISER BOILING VOLUME 

» 

VCLEF=VCIC f*LEF/ LCF 
V3LNL=VJL-L*«.KL/LJL 

* CCMFUTE T H t NONEOIL TIG LENGTH OF THE RISERS 

* 

IDE=LCF* IhF-t-CC)/ ( (HF+ aFF*FFG)-HDD ) 
lwK=LJL* (HF-nJJ) / ( (i-F + aLL*FFGJ-HJJ I 

* 

* CQMOUTE ThE RULING LENGTH DP THE USERS 

* 

LKL=LJL-LJK 
IEF = LCF-LCE 
* 

* 

* 

♦ 

* CD M PUT E THE MASS *?ATE EQUATION FOR STEAM CONDENSING IN 

* TnE CRUM 

MCCNC=56 0 . 93 ^ { om v/ (J mm 4460 .0 ) **U .5-° SAT/ (TS* T +460.U)**0.5)... 

♦ . C2568 

* 

* CC M FMT£ THE RATE EGUATITN FOR DRUM LIQUID MASS 

■* 

comas L=MLL * I 1 . J-XLL )+U F F*( I . C- XFF ) +MCOND+MB30— MCC-MG3 

t 

* COMPUTE THE DRUM LIQUID MASS 



CSMC21I0 
CSMJ2120 
C S"02130 
C SM 0 2 1A0 
CS "02150 
CSM021E0 
CSMC2 1 70 
CS MO 2 180 
C SMC21S0 
CSM02200 
C SM 022 10 
CSMC2220 
CSM02230 
C S "022 AO 
C SM02230 
CS "02260 
CSMC22 70 
C S m C 2 2 80 
CS "C2290 
C SM C2 3C0 
CS M 023 10 
C S V C2320 
CSMC2330 
CS‘-'02JA0 
C S MQ 23 50 
CSMG2360 
CS "023 70 
CSMC23 60 
CS w C23 SO 
CS "02 AQO 
CSMC2A10 
CSM02420 
C SMC2A30 
C S •' t C 2 A A 0 
CS "02A 50 
CS "C2A6C 
CS M 02A 70 
C S"C2A 30 
C SM C2 A SO 
CS "02500 
CS"02510 
C S'M C2 520 
CS "02530 
CSMC25AO 
CS M 02 5 50 
CS "02560 
C SM C 2 5 7 C 
CS v 025e0 
C $"C2 5 90 
C S" C2600 
CS "026 10 
C S "02620 
SS ‘ 02620 
CS " 02 OAO 
CS"C2650 
CSM C2 1 60 
CS -02 6 70 
C S" C2 6EC 
CS" 0 26S0 
CS " 02 7 CO 
CS‘-*C27iO 
CS " 02720 
CS"02 730 
CS" 027AO 
CS "02750 
C S" C2 760 
CS" 02 7 70 
CSM02780 
C S" C2 7 SO 
CS"C2800 
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F CRT RAN PI 



NAVAL p 1ST GR ACL ATE SCHOOL 



F 1 L£ : CSPP 
JhASL s 1NTGRL(DNASLJ .DDPASL) 

*, COMPUTE the PATE EQUATION FCP ENERGY IN the lKU“ liquid 

CJ u OhL="LL* U.J-XLU*HF*aFF*( l.O-XFF ) : *HF +-H CON D^HF G + M3B *HBB . . . 

— FCC*HCC-*v>vj*hGG 

* CO-NP'JTE THE ENERGY IN THE DRUM LIQUID 

D^CHL= INTG C L ( QHQHL C » OOMDHL ) ' 

* COMPUTE THE ENTHALPY CF THE DRUM LIQUID 

* 

OH=OMDHL /CMASL 
* 

* 

* CUMPITE THE DRUM SPECIFIC VOLUME 

* 

CRYUNO=PCU* iFLl+ v FF) 

PROCEDLRE R I SE 2 =F l lTRA ( C RYUND) 

R IS61=CtL i Y ( ^ 50 »R 1ST I M , CR Y LND ) 

R I SE2 =CRYU AD 

IF ( VALVE. GT. . 5 1 ) GO TO 55 
R IS E= R IS E2 
GC TC 57 

55 ■ R I S t = R IS E l 

57 CCM1ME 

EN C°POC ECU’ fc 

CCMOV= ( ( v F F+ V L L~N 69} »VF*’C I *( M FF +mlL ) *VV *MB9 0*VB8-( MCC*-MGG) *VF <-MD ONO . . 

* VPG-«I SE * V V ) 

CM0V= I.NTGRL ( DND'/O *DDMDV ) 

* 

* COMPETE THE DRU > LEVEL 

* 

LEVE L=(Dmcv-VCLoR' v /2.0 }/ ( IRT ^ G >1 * CS t M D'-< ) 

* COMPETE THc DOWNCOMER ENTHALPY 

* 

HCD= I ( FFF + F’LL-HB5)~HF-*-VE' , < : FB3}/ t FFF + ^LL) 

FGF= FC D 
FCC=FCC 
FCC=HCD 
FGC-= FGH 
* 

* COMPUTE THE CCwr^CCMER SPECIFIC VOLUME A NO DENSITY 

* 

VC 0=( ( FFF ♦FLL-FBB) * VF + FEB* VBB) / ( FFF+*LL‘ 

®HOCD= l.O/VCO 
RFCGh = RH CCD 
* 

* fO‘ ,D UTE thE SATURATION PRESSURE AND TEMPERATURE 

* CCPPtS PCNCINl TC T Ft CDWNCRMER ENT H AL P Y 

* 

PS AT=EXP ( ( ALOC-IFSAT !-<♦ .Aj7C3 )/.2oA 5 2J - 
T SAT = E >Pt {.22I5 L*ALO^(PSA t )+4./7123 ) ) 

* 

* CCMF’ITf Tht ElTFALPY CF TFE LIQUID ENT ER I NG TFE 

* MAIN BANK RI SER 

4 

H J J=HGH+ C9 /PHH 

ft 

* CC v PUT £ THE RATE EQUATION FDR THE MAIN BANK AND 

* SCREEN RISER fUEE METAL TEMPERATURES 

< 

OTVVMCI-C?)/ ( NAS$V*CPV] 

OTYV=(Q5-QO /( MASSY*CPY) 

* CCMPLTE THE SCREEN AND M A 1 N BANK RISER TUBE 

« METAL tfmpef. AFURES 

* 



CSMC28IO 
C S v 023 20 
CS MJ 2 8 30 
CS V C2840 
CS M C2 8 50 
CS ‘-'C 23 60 
CSMC28 70 
CS M 02 3 SO 
C5MC2890 
CSMC2SC0 
C 3 M 029 10 
C SM C2920 
CS M C 29 20 
CSMC2940 
C S M C 29 50 
CS v 029 60 
CSMC29 70 
C SM 0 29 80 
CS M 029 90 
C S v C 2 C CO 
C S • G 3 0 10 
CS MC 3020 
C S M C A C30 
CS MO 3040 
CSMQ3050 
CSMC2C6C 
CS M G3 0 70 
C SM C3C80 
CS v 0 3090 
C 3 *0 3100 
C 5 “ C 2 1 10 
C S v 0 3 1 20 
C S^C3130 
CS m o 3 1 40 
CS V C3 1 50 
C SMC J 160 
CS V 03170 
CS*C3180 
: S N< 0 3 1 9 0 
CS *03200 
CSMC2210 
C 3 ‘AC 22 20 
CS MO 32 30 
C SM C 32 40 
CS M 0 2 2 50 
C Sv 032 60 
C 3 '•* C 3 2 70 
CS VO 3 2 80 
C SMC 32 90 
CS* 0 3500 
CS * Cj 3 10 
C SMC 2 3 20 
CS* Jj>330 
C S VC 3 3 40 
: S'* G 33 50 
CS Vo J360 
C SM C 2 2 70 
CS *•' 033 80 
CS M 3 3 390 

C 5* C34C0 
C3 V C3410 
C SMC24 20 
CSM C 3420 
C S v C 34 40 
C SMC3450 
CSM 03460 
C S M C34 70 
C S* C2480 
CS *03490 
CS VC35CC 
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f-ILfis CS^P 



FORTRAN PI 



NAVAL POSTGRADUATE SCHCOI 



rW=INTGRL( TVVC,DTVV) 
r rr= 1NTGPL (TYYO, CTYYi 

* COMPUTE THE HEAT TRANSFER FOR ThE DESUPERHEATER 

* -STEA M TO TUBE FETAL 

* 

09 =C N p *MNN*( TNN-TPP J 



* 

* 

♦ 

* 



CONFUTE THE TEMPERATURE OF THE STEAF LEAVING THE 
CESUPEPHEATER 



TPP= ( TNN-T 11 ) / ( E XP { KNP / ( CN° *MNN**0*2 ) ) ) + TZZ 



SET THE DOWNCOMER ENTRANCE AND EXIT T EM P ERATIJR ES 
EQUAL TO THE SATURATION TEMPERATURE C ORRE SPONDI NG 
TC CRUF ENTMLPY 



COMPUTE THE ENERGY TRANSFER FOP THE DESUPERHE AT ER 
TUBE M ET AL TC DRUM L I OU I D 



COMPUTE THE DESUPERHEATER TUBE *ETAL TEMPERATURE 



♦ 

* 

* 

* 

* 

TGG=TS *T 
7HH=TSAT 
TC C= TS AT 
TDC=TS AT 
* 

♦ 

♦ 

* 

Q9 9=KZ*( TZZ-Thh) 

♦ 

* COMPUTE THE RATE EQUATION FOR D ESU PERHt AT ER TIJEE METAL 

4 TEMPERATURE 

* 

CT ZZ= ( Q9-Q 99 ) / ( MAS$Z*CPZ) 

* 

* 

4 

TZZ = INTGRU T ZZ0tCT2Z) 

♦ 

♦ COMPUTE ThE ECONOMIZER ENERGY TRANSFER-TUBE METAL 

♦ TO FEED WATER 

4 

CS = MAA*CA E* (T EE-TAA ) * 

♦ 

* compltf the feed temperature at cutlet 

* CF ECONO M I 2 ER 

* 

T BB= (TAA-TXX) / ( EXP ( KX/ ( C A ° *MA £ ** C. 2) ) ) +TXX 

4 

♦ rOMPL'I’E THE RATE EQUATION FOR THE ECONOMIZER TUEE 

♦ TEMPERA T URE 

4 

OT XX = ( 08-07)/ ( M A $ 3 X *C P X ) 

♦ 

* COMMUTE THE ECONOMIZE R TUEE M £T AL TEMPERATURE 

* 

T>>=I NTUPL (TXXC f CTXX J 
* 

* CORFU'S THE SJPERHEATEC STEAM OUTLET TcMPERATLRE 

* 

TNN= (TMM-TWW i/( EXP(RW/(CMN *MNN**C. 2) )) +TWW 
* 

* COMPUTE THE SUPERHEATER ENERGY TRANSFER-TUBE w £TAl 

* TO STEAM 
C4=C FN*MFM* i TNN-TMM) 

4 

* C0 M pu t E THE RATE EQUATION FOR SUPERHEATER TUBE METAL 

* TEMPERATURE 

4 

CTWW = (Ci-04 )/ ( N AS S W*C PW ) 

4 

4 COMPUTE THE SU° ERH EA TER TUBE METAL TEMPERATURE 



CSM03510 
CS MO 35 20 
CSMC3520 
CS *03540 
C SMC 25 50 
CSMC3560 
CS M 035 70 
C SMC35EC 
CSMC3590 
C S M C36 00 
CSMC36I0 
C S *036 20 
C S MO 3 63 C 
CS M 03640 
CS M 036 50 
C SM C 36 60 
CSMC3670 
CSMO260O 
C SM C3690 
CS M 03 7 00 
C S MQ3 7 10 
CSM02720 
CS m J3 7 30 
C SMC3740 
C S M 03 7 50 
CS M03760 
C SM C 3 7 70 
CSM037EO 
CSM03790 
C SMC3EC0 
CS M038I0 
C S M C 2 8 20 
C S C 3 6 20 
CS M 03 840 
C SMC3650 
CS M 03 860 
C SMC3870 
CSMC3860 
CS MO 3 090 
C S 39 00 
CSM039 10 
CS *0 39 20 
CSMC3930 
C S m C 3 9 4 0 
CS '*03950 
CSMC2960 
CS *0 29 70 
C SM C 39 80 
CSMC2S90 
CS *04000 
C S M C4C10 
CS“04020 
CSM04030 
C SM C 40 40 
CS MO 40 50 
CSMC4C60 
CS V 04 0 70 
CS M C-+ G 80 
C S M C 4 0 90 
CS s * 04 1 00 
CS*C4I 10 
C S* Oh 120 
CS MO 41 30 
C S * C4> 1 40 
C SM C 4 1 50 
CSM04160 
CS M C4I 7C 
CS* 04 1 80 
CS MC4I90 
CSMC42C0 
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HIE: CS vp 



T *i,= I NTGRL< TV* k>C tOTWVi) 

» ccpfute the implicit equation 

• FlC* SATE FFCR E01LER 



FOR TOTAL 



M»F = I FFLI*VN 3, NMC’E, VMM IFP) 

C n R PUT E THE SJPFRhEATEP OUTLET PPESSU 3 E 

oHN=SO p T ( OS AT * + l . 0 - 2 . 0 *K 0 N 2 *K 0 N 1 *mmm** 2 . 0 ) 

II = p r'MKCN4* VALVE 

VMM£ I = TJ 0 

NNH kplMf/M J ] + M MM III 



STEam 



* FOLATE the flue gas mass flow rates 

* 

4RR=M0Q 

•*‘SS=MRR 

► T T= MS S 
MLL=MTT 

4 

* SET RISER INLET CEN'S ITY EQUAL TO 

* SATURATED LIQUID DENSITY 

* 

PH CEE = RH CF 
R HOJ J= RH CF 
C HCCC= SHOE E 
RH CK N = RH C J J 

* . 

* ECU i T E COWNCCMfcR FL r W RATES TO 0 IS E p EN T 3 ANO E FLOW 

* RATES ANO SET RISER ENTRANCE FLCW PAT£S r j THE 

* FLOW RATES AT THE INITIAL TIME 

* 

VHH=MGG 
>CC= V DC 
PGG= NJ J 
MDD=MFFC 
v J J= v llo 
« 



* COMPUTE THE DERIVATIVE JF AVERAGE RISER DENSITY 

* 

* 

* -COMPUTE THE STEAM MASS RATE EQUATION FOR THt STEAM 

* DRUM 

DDSTM= AFF*mFF + >LL* m LL- VCOND-mmm 

* 



* COMPUTE ThE CRUM STE A v MASS 

4 



CSTM= INTER L ( CS T v 0 t DOST M J 
* 



* COMPUTE THE VOLUME OF STEAM IN THE STEAM DRUM 

* 



VnLSTM=VQLCRM-DMA SL/RHOF 
* 

* COMPUTE THE DENSITY OF STEAM IN THE STEAM C P UM 

* 

RHCSTM = C'THVCLST-M 

rhomm=rhcstm 

V'/mm = 1 .0/RHOMM 

vfgmm=vvem-vf 

* 

* COMPUTE THE STEAM DRUM STEAM OUTLET PRESSURE 



p°ESSM=524.(/< VFGMM+. 1) 

PROCEDURE PMM = F1LTP5 (FPESSM) 

IF(P°FSSm.LT. 0. CJCALL DE3UG (3,0.0) 

FVv=pfe$$v 

ENDPRCCkCLPE 

* COMPUTE THE STEAM DRJM STEAM OUTLET TEMPERATURE 



CSM C42 10 
CS MC4220 
C S M C42 30 
CS M Cw240 
CSM04250 
C SMC<*2tO 
CS v 042 70 
C S M C 4 2 8 0 
CS M C42S0 
CS M 04300 
C SM C4310 
CSM 04320 
C SM04330 
CSM 04340 
CS M 043 50 
C SMC4360 
CSMQ4370 
CS MJ 43 60 
C SMC 42 SO 
CSMC4400 
CS MQ4410 
C SMC4420 
CS’ 1 04430 
CSNC4440 
C SM C 44 50 
CSMQ4460 
CSM04470 
C SM C4480 
CS M044S0 
C S C 4 5 G 0 
CSMG45 10 
CSM04520 
CSM C4530 
CSM 04 540 
C S M 04550 
C SMC 4 560 
CSM04570 
CSM04590 
C SM C45S0 
CS “04600 
C S M C 4 1 10 
CSMC4620 
CS m C46 30 
C SM C464C 
CS v 04650 
CSMC4b60 
CSMC4670 
CS M046 80 
CS VC4 6S0 
CS m 04 7 00 
C S“G4710 
CSM C4 7 20 
CS “04730 
CSM04740 
C S M C47 50 
CS M 04 760 
CSMC4770 
CSMC4760 
CS “04790 
C SM C4 8C0 
CS MC4810 
C SMC4820 
C SM04830 
C S M 04S 40 
C S“C4650 
CSM 04 860 
C SMC4870 
CSMC 4880 
CS “04690 
CSMC4900 
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ML£: csmP 



FORTRAN PI 



NAVAL PQSTGRAQLATE SCHCOL 



rM«=cX D ( . 22 IS 1*AL0G(PMM) +4.77123) 

* SJL \E FOR THE ECCNOMI 2ER FEED OUTLET ENTHALPY 

HSE« (>'eB*HAA0 + C8)/MBe 

* fEECRATE EQUATION 

* 

*BB=M80O 

* 

* two phase flow *ultiplier equations 

* 



PG P AV E=2 4 . 7°4 *X F F * *2 .0— 6.50t6*XFF+.9776 
RGPA VK=24. 794* XL L** 2.0- t. 57 6 6* XL L+ . 9 776 



- A CL £= ii 


j.45fc-* 


x xFF **2.0 


+ lo. 4944*XFF— , 


, C0007 




RACL J=l 5.456-+ 


*XLL **2 . J 


+ 13 .4 94 “♦ *< L L — , 


.00307 




5FRI CE=- 


■ 34 . Ci 


22* Xr F ** 2 


. 0+23. 7164* XFF +.8734 




FFFi CK=- 

X 


.0 a 


2 2 *XLL **4 


.C+2 j. 7164*XLL ♦. 8734 




* 

* 


STATE POINT 


EOLATIONS 






w 

nS^T^OH 












br= HSA T 












b F 1= 9 2 2 « 


.15-0 . 


435 16 ^PS A 


T + 1.717E-0 4*°SAT**2.C- 


4.219E-0 


*PSAT**3 


• 0 








FV = HS at * 


M- FG 




' 






FrtCF =63 . 


► 8-0*0 


L731* T S AT 


+1 .132E-0 5 *TSAT**2 .0-6 


. 786E-03 



C S M C 49 10 
CS ^049 20 
CS*C4930 
CS*'* C49 40 
CS *04950 
CSMC4560- 
CS* C4970 
CSR049RO- 
C S M C4550 
CS *'0 5000 
CSHC5010 
C S M C 5C20 - 
CS *05030 
CS *C5C40 
CS**a50 50 
CS ^05060 - 
C SMC5C70 
CS '<05080 - 
CSF0509a - 
C SMC 5 100 
CS *05110 - 
C S M 35 1 20 
CSMC5I30: 
CS»'05140 
CSFC5I50* 
CS *'05 160 
CS *05170 
CSMC51E0 
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Gv Anticipated Performance 



APPENDIX C 



a of Operation * Per Cent 
i o l Steam Generated lb/hr 
qrheated Steam Ibs/hr. 
uC'perheated Steam Ibs/hr. 

Bo r Drum Pressure psig 
ujrhcater Outlet Pressure psig 
Sujrheater Outlet Temperature °F 
De:perheated Steam Outlet Pressure psig 
De.'pcrheated Steam Outlet Temperature °F 
Eaomizer Inlet Pressure psig 
Eaomizer Inlet Temperature °F 
Eaomizer Outlet Pressure psig 
xomizer Outlet Temperature °F 
Castg Air Inlet Temperature °F 
Toil Air Flow Ib/hr 
Toll Oil Flow ib/hr 
Aji ;ipated Efficiency % 
uU;anteed Efficiency % 

Raotion and Unaccounted for Losses % 

Extss Air % 

Cann Dioxide % 

Naiber of Burners in Operation 
Ihittle or Non-Throttle of Air Doors 

Dra Loss - Total Inches Water 
Ulugh Double Casing 
rhugh Burner Register 
Fhiagh Boiler Sc Superheater 
ujgh Economizer 

Jas’emperature Leaving Superheater Screen °F 
Jas'emperature Leaving Superheater °F 
^as'emperature Leaving Main Bank °F * 
Jas’emperature Leaving Economizer °F 
eaRelease KB/Hi/Sq. Ft. Radiant Heat Absorbing 
>ur ce 

03 Release KB/Hr/Sq. Ft. Total Heating Surface 
ieaRelease KB/Hr/Cu. Ft. Furnace Volume 
urice Ileal Absorption KB/1 Ir ./Sc? . Ft 
ea Absorption First Water Scieen Row KB/Hr/Sq. 
Ft. (Max.) 

lea Absorption Maximum KB/llr./Sq. Ft 
(Furnace Screen) 



Endurance 


Design 

Rated 

Full 


Maximum 

Intermittent 


Endurance Boiler 


2 Boilers 


Power 


Power 


1 Boiler 


Overload 


56 


83 


97 


1 14 


120 


195,000 


290,000 


340,000 


400.000 


420.000 


155,000 


235,000 


265.000 


300,000 


320.000 


40,000 


55,000 


75,000 


100,000 


100.000 


690 


690 


690 


690 


690 


665 


b30 


610 


580 


570 


912 


914 


980 


898 


895 


651 


606 


558 


480 


470 


635 


659 


680 


690 


690 


715 


726 


733 


749 


754 


250 


250 


250 


250 


250 


703 


705 


703 


708 


710 


355 


372 


372 


377 


380 


100 


100 


100 


100 


100 


242,229 


368.092 


433,589 


513,154 


540,952 


14,333 


21,781 


25,656 


30,364 


32,009 


85.7 


84.4 


83.8 


83.0 


82.7 


85.7 


84.4 


83.8 


83.0 


82.7 


.99 


.99 


.95 


1.01 


.98 


15.0 


15.0 


15.0 


15.0 


15.0 


13.0 


13.0 


13.0 


13.0 


13.0 


6 


6 


6 


6 


6 


N.T. 


N.T. 


N.T. 


N.T. 


N.T. 


13.74 


33.98* 


47.96 


69.51 


76.22** 93.48 


1.33 


3.09 


4.28 


6.00 


6.67 10.22 


2.60 


6.00 


8.50 


1 1.50 


13.00 15.88 


4.83 


12.29 


18.14 


25.61 


28.15 38.32 


4.98 


12.60 


17.04 


26.40 


28.40 29.06 


2471 


2594 


2642 


2686 


2699 


1826 


1961 


202 J 


2085 


2103 


694 


774 


812 


856 


872 


373 


427 


453 


484 


498 


447 


680 


800.6 


948 


999 


15.5 


23.6 


28.0 


32.9 


34.7 


183.9 


279.5 


329.2 


389.6 


410.S 


128.4 


171.0 


191.1 


213.19 


221.5 


210 


231 


255.3 


28S 


256 


210 


231 


255.3 


288 


256 



r; Losses (Full Power) Based On 24 1 Cu. Ft. of 100°F Air/Lb. of Fuel Oil 
*D ft Losses (Overload) Based On 241 Cu. Ft. of 10C°F Air/Lb. of Fuel Oil 
aft Losses (Overload) Based on 260 Cu. Ft. of 68°F Air/Lb. of Fuel Oil 
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II. Anticipated Metal Temperature Degree F. 





Endurance 


Design 

Rated 


Maximum 

Intermittent 


Endurance 


Boilei 




2 Boilers 


Full Power 


Power 


1 Boiler 


Overload 


Water Screen Tubes Outside 


644 


656 


661 


667 


680 


Water Screen Tubes Inside 


554 


554 


554 


554 


555 


Superheater Tubes Outside - Maximum 


1022 


1046 


1049 


1048 


1048 


Superheater Tubes Inside - Maximum 


1044 


1020 


1019 


1014 


1012 


Maximum Inner Casing Temperature 


SOO 


825 


835 


850 


870 


Maximum Outer Casing Temperature 
at way of Structural ties 


350 


350 


350 


350 


350 


Maximum Outer Casing Temperature 


145 


145 


145 


145 


145 



Tube Data 


O.D. 


M.W.T. 


No. 


Material 


Side Wall and Roof 


2” 


.134” 


71 


MlL-T-16286 CL. A 


Rear Wall 


2” 


.134” 


54 


MIL-T-16286 CL. A 


Front Wall 


2” 


.134” 


22 


MIL-T-16286 CL. A 


Screen Bank 


2” 


.134” 


102 


MIL-T-16286 CL. A 


Superheater 


1.5” 


.120 


22 


MIL-T-16286 CL. E 


Main Bank 


2” 


.134” 


34 


MIL-T-16286 CL. ;v 


Main Bank 


1” 


.085” 


2808 


MIL-T-16286 CL. A 


Economizer 


2” 


.180” 


182 


M1L-T-16286CL. a 


Desuperheater (in Water Drum) 


2” 


.220” 


6 


SA-268 TP430 


Downcomers . 

* 


8 5/8” 


.483” 


1 


Schedule 80 Pipe 
ASME SA-106-B 


Downcomers 


10 3/4” 


.519” 


5 


Schedule 80 Pipe 
ASME SA-106-B 


Downcomers 


12 3/4” 


.601 


2 


Schedule 80 Pipe 
ASME SA-10M1 


Risers 


6” 


.500” 


7 


ASME S A- 106-11 



““These tubes may be MlL-T-16286, Class A (Seamless) or MILT-17188 (Seamed) electric resistance welded. 
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